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Abstract: We present a detailed study of the finite-temperature behaviour of the LARGE 
Volume type IIB flux compactifications. We show that certain moduli can thermalise at 
high temperatures. Despite that, their contribution to the finite-temperature effective 
potential is always negligible and the latter has a runaway behaviour. We compute the 
maximal temperature T max , above which the internal space decompactifies, as well as 
the temperature T* , that is reached after the decay of the heaviest moduli. The natural 
constraint T* < T max implies a lower bound on the allowed values of the internal volume 
V. We find that this restriction rules out a significant range of values corresponding to 
smaller volumes of the order V ~ 10 4 /!?, which lead to standard GUT theories. Instead, 
the bound favours values of the order V ~ 10 15 Zf , which lead to TeV scale SUSY desirable 
for solving the hierarchy problem. Moreover, our result favours low-energy inflationary 
scenarios with density perturbations generated by a field, which is not the inflaton. In 
such a scenario, one could achieve both inflation and TeV-scale SUSY, although gravity 
waves would not be observable. Finally, we pose a two-fold challenge for the solution of the 
cosmological moduli problem. First, we show that the heavy moduli decay before they can 
begin to dominate the energy density of the Universe. Hence they are not able to dilute 
any unwanted relics. And second, we argue that, in order to obtain thermal inflation in 
the closed string moduli sector, one needs to go beyond the present EFT description. 
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1. Introduction 

The low energy effective action of string compactifications on Calabi-Yau (CY) 3-folds 
typically has a large number of uncharged massless scalar fields with a flat potential, called 
moduli. This has been a long-standing problem for string phenomenology, as the values of 
those moduli determine the various parameters (like masses and coupling constants) of the 
four-dimensional effective description. Hence, the presence of these massless scalars with 
effective gravitational coupling would lead to unobserved long range fifth forces, as well as 
a lack of predictability of the theory. 

The last decade has seen a lot of progress towards the resolution of this problem. A 
major ingredient in these developments was the realization of ||, |2|, || that nonzero back- 
ground fluxes induce potentials for some of the moduli. In fact, in type IIA all geometric 
moduli can be stabilized in this way j3j. In type IIB, on the other hand, one also needs 
to take into account various perturbative and non-perturbative effects |5|, ^] . Hence it may 
seem that IIA compactifications are under better control. The reason this is not so is that 
the backreaction of the fluxes on the geometry is more severe in type IIA than in type IIB. 
As a result, generically in type IIA one has to consider as internal spaces manifolds with 
SU(3) x SU(3) structure (see 0] for a comprehensive review). The latter are mathemati- 
cally much more involved than a Calabi-Yau and so are, in principle, much harder to study. 
In contrast, in type IIB there is a huge class of solutions, such that the backreaction of the 
fluxes is entirely encoded by a warp factor. Naturally then, it is within this class of IIB 
CY flux compactifications that moduli stabilization is best understood at present. 

An excellent example of IIB compactifications with stabilized moduli, which is very 
appealing both for particle physics phenomenology and for cosmology, is given by the 
LARGE Volume Scenarios (LVS) originally proposed in §. According to the general 
analysis of |8|, in these compactifications, a' and g s corrections are combined with non- 
perturbative effects to generate a potential for the Kahler moduli, whereas the background 
fluxes induce a potential for the dilaton and the complex structure moduli. Unlike the 
KKLT set-up ||, the moduli stabilisation is performed without fine tuning of the values 
of the internal fluxes and the CY volume is fixed at an exponentially large value (in string 
units). As a consequence, one has a very reliable four-dimensional effective description, as 
well as a tool for the generation of phenomenologically desirable hierarchies. 

The exponentially large volume minimum of LVS is AdS with broken SUSY, even before 
any uplifting. In contrast, in KKLT constructions the AdS minimum is supersymmetric 
and the uplifting term is the source of SUSY breaking. In both cases, however, in addition 
to this minimum there is always the supersymmetric one at infinity. The two minima are 
separated by a potential barrier V^, whose order of magnitude is very well approximated 
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Figure 1: The effective potentiai V versus the voiume modulus <fi for a typical potential of KKLT 
or LARGE Volume compactifications. The different curves show the effect of various sources of 
energy that, if higher than the barrier of the potential, can lead to a decompactification of the 
internal space. 



by the value of the potential at the AdS vacuum before uplifting. As is well-known, the 
modulus, related to the overall volume of the Calabi-Yau, couples to any possible source of 
energy, due to the Weyl rescaling of the metric needed to obtain a 4D supergravity effective 
action in the Einstein frame. Thus, in the presence of any source of energy, greater then the 
height of the potential barrier, the system will be driven to a dangerous decompactification 
limit. For example, during inflation the energy of the inflaton ip could give an additional 
uplifting term of the form AV(ip,V) = V(ip)/V n for n > 0, that could cause a run-away 
to infinity 11. Another source of danger of decompactification is the following. After 
inflation, the inflaton decays to radiation and, as a result, a high-temperature thermal 
plasma is formed. This gives rise to temperature-dependent corrections to the moduli 
potential, which could again destabilise the moduli and drive them to infinity, if the finite- 
temperature potential has a run-away behaviour. The decompactification temperature, at 
which the finite-temperature contribution starts dominating over the T = potential, is 
very well approximated by T max ~ V^ 4 since Up ~ T 4 . Clearly, T max sets also an upper 
bound on the reheating temperature after inflation. The discussion of this paragraph is 
schematically illustrated on Figure [l|. 

On the other hand, if, instead of having a run-away behaviour, the finite-temperature 
potential develops new minima, then there could be various phase transitions, which might 
have played an important role in the early Universe and could have observable signatures 
today. The presence of minima at high T could also have implications regarding the 
question how natural it is for the Universe to be in a metastable state at T = 0. More 
precisely, recent studies of various toy models 11, [T^, 13, Q have shown that, despite 



the presence of a supersymmetric global minimum, it is thermodynamically preferable for 
a system starting in a high T minimum to end up at low temperatures in a (long-lived) 
local metastable minimum with broken supersymmetry. Similar arguments, if applicable 
for more realistic systems, could be of great conceptual value given the present accelerated 
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expansion of the Universe. 

For cosmological reasons then, it is of great importance to understand the full structure 
of the finite temperature effective potential. We investigate this problem in great detail 
for the type IIB LVS of Q. Contrary to the traditional thought that moduli cannot 
thermalize due to their Planck-suppressed couplings to ordinary matter and radiation, we 
show that in LVS some of the moduli can be in thermal equilibrium with MSSM particles 
for temperatures well below the Planck scale. The main reason is the presence of an 
additional large scale in this context, namely the exponentially large CY volume, which 
enters the various couplings and thus affects the relevant interaction rates. The unexpected 
result, that some moduli can thermalize, in principle opens up the possibility that the 
finite temperature potential could develop new minima instead of just having a run-away 
behaviour as, for example, in fll5| . However, we show that this is not the case since, for 
temperatures below the Kaluza-Klein scale, the T-dependent potential still has a run-away 
behaviour. Although it is impossible to find exactly the decompactification temperature 
T m ax, as it is determined by a transcendental equation, we are able to extract a rather 
precise analytic estimate for it. As expected, we find that T max is controlled by the SUSY 
breaking scale: T^ ax ~ m^ 2 Mp. This expression gives also an upper bound on the 
temperature in the early Universe. We show that this constraint can be translated into a 
lower bound on the value of the CY volume, by computing the temperature of the Universe 
T*, just after the heaviest moduli of LVS decay, and then imposing T* < T max . 

Our lower bound implies that, for cosmological reasons, larger values of the volume 
of the order V ~ lO 15 /^, which naturally lead to TeV scale super symmetry, are favoured 
over smaller values of the order V ~ 10 4 ^, which lead to standard GUT theories. More 
precisely, what we mean by this is the following. Upon writing the volume as V ~ 10 x 
and encoding the fluxes and the Calabi-Yau topology in the definition of a parameter c, 
we are able to rule out a significant portion of the (x, c)-parameter space that corresponds 
to small x (for example, for c = 1 we obtain x > 6). This is rather intriguing, given that 
other cosmological considerations seem to favour smaller values of the volume. Indeed, 
the recent inflationary model of [16] requires V ~ 10 4 /g, in order to generate the right 
amount of density perturbations. Despite that, our lower bound on V does not represent 
an unsurmountable obstacle for the realization of inflation. The reason is that the Fibre 
Inflation model can give rise to inflation even for large values of the volume. Hence, a 
modification of it, such that the density fluctuations are generated by a curvaton-like field 
different from the inflaton, would be a perfectly viable model with large V. The large value 
of V would imply a low-energy inflationary scale, and so, in turn, gravity waves would not be 
observable. However, it is likely that both inflation and TeV scale SUSY could be achieved 
at the same time, with also the generation of a relevant amount of non-gaussianities in the 
CMB, which is a typical feature of curvaton models. 

On the other hand, we pose a challenge for the solution of the cosmological moduli 
problem, that the overall breathing mode of LVS with V > 10 10 /f is afflicted by [ IT | . This 
is so, because we show that unwanted relics cannot be diluted by the entropy released by 
the decay of the heaviest moduli of LVS, nor by a low-energy period of thermal inflation. 
More precisely, we show that the heaviest moduli of LVS decay before they can begin to 
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dominate the energy density of the Universe and, also, that in order to study thermal 
inflation in the closed string moduli sector, it is necessary to go beyond our low energy 
EFT description. 

The present paper is organized as follows. In Section ||, we review the basic features 
of the type IIB LVS. In Section ||, we recall the general form of the effective potential at 
finite temperature and discuss in detail the issue of thermal equilibrium in an expanding 
Universe. In Section ||, we derive the masses and the couplings to visible sector particles 
of the moduli and modulini in LVS. Using these results, in Section || we investigate moduli 
thermalization and show that, generically, the moduli corresponding to the small cycles 
can be in thermal equilibrium with MSSM particles, due to their interaction with the 
gauge bosons. In Section ^, we study the finite temperature effective potential in LVS. 
We show that it has a runaway behaviour and find the decompactification temperature 
Tmax- Furthermore, we establish a lower bound on the CY volume, which follows from the 
constraint that the temperature of the Universe just after the small moduli decay should 
not exceed T max . In Section we discuss some open issues, among which the question why 
thermal inflation does not occur within our approximations. Finally, after a summary of 
our results in Section ^, Appendix [A] contains technical details on the computation of the 
moduli couplings. 

2. Large Volume Scenarios 

The distinguishing feature of the type IIB Large Volume Scenarios (LVS) of || is that, 
in addition to the non-perturbative effects considered in ||, they also take into account 
a' corrections, which lead to moduli stabilisation at an exponentially large volume of the 
internal manifold. In Q, the Calabi-Yau 3-fold was assumed to have a characteristic 
topology with one exponentially large cycle and several small del Pezzo 4-cycles. Including 
string loop corrections, as in ||, extends these scenarios to a larger class of Calabi-Yau 
manifolds which can also have fibration structure. The exponentially large volume allows 
to explain many hierarchies observed in nature and guarantees that the low-energy effective 
field theory is under good control. 

In this Section we summarise the basic ingredients of LVS. We begin by recalling 
necessary material about type IIB flux compactifications. After that we turn to the relevant 
perturbative and non-perturbative effects and the resulting scalar potential with a LARGE 
volume minimum. 

2.1 Type IIB with fluxes 

We will be interested in type IIB CY orientifold compactifications with background fluxes, 
which preserve Af = 1 supersymmetry in 4D ||. The effective low-energy description 
is then given by a d = 4, J\f = 1 supergravity characterised by a Kahler potential K, 
a superpotential W and a gauge kinetic function f a b, where the indices a, b run over the 
various vector multiplets. In particular, the scalar potential of this theory has the standard 



- 5 - 



form: 

|2- 



V = e R / M P ^DiWD^W - 3^r) , (2.1) 



where i,j run over all moduli of the compactification. Generically, the latter consist of 
the axio-dilaton S = + iCq, Kahler (T-moduli) and /i2,i complex structure moduli 
([/-moduli). The tree level Kahler potential has the following form: 



K, 



tree 



M 2 p 



In (S + S) -21nV-ln j -i J 



n Aft , (2.2) 



CY 



where V is the Einstein frame CY volume, in units of the string length l s = 27rv / c7, 
and ft is the CY holomorphic (3,0)-form. Note that the the T-moduli enter Kt ree only 
through V and the [/-moduli only through ft. For later purposes, it will be useful to 
recall a couple of relations regarding the Kahler moduli. Expanding the Kahler form 
in a basis {Di} of H 1,1 (CY,7<) and considering orientifold projections such 
that h^i = =>■ = hi i, we obtain: 

V = ll J^J^J=lk ij kt i t j t k , (2.3) 
6 .Icy 6 



where are related to the triple intersection numbers of the CY and the t 1 are 2-cycle 
volumes. The volumes of the Poincare dual 4-cycles are given by: 

Ti = % = \ / D '< AJAJ = \ k i3kt 3 t k . (2.4) 

CY 

Finally, the scalar components of the chiral superfields, corresponding to the Kahler moduli, 
that enter the 4D effective action are T{ = n + ibi, where the axions bi are the components 
of the RR 4- form C4 along the 4-cycle Poincare dual to Di . Obviously, from ( [2.3| ) and (p.4| ) 
one can express V in ( |2.2| ) as a function of Tj = \{Ti + Tj). 

Now, turning on background fluxes G3 = -F3 + 1SH3, where F% and H3 are respec- 
tively the RR and NSNS 3-form fluxes of type IIB supergravity (for recent reviews on flux 
compactifications, see p8|), generates an effective superpotential of the form |l( : 



Wtree ~ / G3 A fi. (2.5) 

CY 

As a result, one can stabilise the axio-dilaton S and the [/-moduli. However, the Kahler 
moduli Tj do not enter Wtree an d therefore remain massless at leading semiclassical level. 
One can obtain an effective description for these fields by integrating out S and U. 2 Then 
the superpotential is constant, W = (Wtree) = Wq, and the Kahler potential reads K = 
K cs - In (2 /g s )+K Q with 

K = -2lnV and e~ Kcs = l-i J QA^ . (2.6) 



1 We neglect the effect of warping, generated by non-zero background fluxes, since we will be considering 
CY compactifications with large internal volume. 

2 For more details on the consistent supersymmetric implementation of this procedure, see E^|. 
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Before concluding this Subsection, let us make a couple of remarks. First, note that the 
background flux G3 may or may not break the remaining 4D ]\[ = 1 supersymmetry, 
depending on whether or not D a W = d a W + Wd a K vanishes at the minimum of the 
resulting scalar potential. Second, the Kahler potential Kq satisfies the no-scale property 
Kq diKod^Ko = 3, which implies that the scalar potential for the Kahler moduli vanishes 
in accord with the statement above that the fluxes do not stabilise those moduli. 



2.2 Leading order corrections 

As we recalled in the previous Subsection, at leading semiclassical order the scalar potential 
for the T-moduli vanishes. So, unlike the situation for the S and [/-moduli, in order to 
stabilise Tj one has to consider the leading order corrections to the tree level action. The 
first such corrections to be studied were non-perturbative contributions to W. Recall that 
there is a non-renormalisation theorem forbidding W to be corrected perturbatively. On 
the other hand, the Kahler potential K does receive perturbative corrections both in a' 
and in g s . Therefore, non-perturbative effects are subleading in K and we shall neglect 
them in the following. Let us now review briefly all these kinds of corrections. 



2.2.1 Non-perturbative corrections 

Non-perturbative corrections to the superpotential can be due to Euclidean D3 brane (ED3) 
instantons wrapping 4-cycles in the extra dimensions, or to gaugino condensation in the 
supersymmetric gauge theories located on D7 branes that also wrap internal 4-cycles. The 
superpotential that both kinds of effects generate is of the form: 3 

W = ^t(w + ^2^e- a A, (2.7) 



47T 

\ 't / 

where the coefficients correspond to threshold effects and, in principle, can depend on 
U and D3 position moduli, but not on Tj. The constants a, are given by dj = 2n for 
ED3 branes, and by a% = 2ir/N for gaugino condensation in an SU(N) gauge theory. 
Note that one can neglect multi-instanton effects in (^]7]), as long as Tj are stabilised such 



that OjTj 3> 1. From (24), the above superpotential leads to the following Tj-dependent 
contribution to the scalar potential (up to a numerical prefactor and powers of g s and Mp, 
that we will be more precise about below): 



SV{n P ) = e K °K 



0,-Af aiA ie -( a ^ +aiT 



1^3 

djAje-^WdiKo + aiAie- aiTi WdjK ) . (2.8) 



3 The prefactor in (2.7|) is due to careful dimensional reduction, as can be seen in Appendix A of [G9[ 



However, the authors of [M define the Einstein metric via g^v,, 



,W-M)/2 



so that it coincides with 



the string frame metric in the physical vacuum. On the contrary, we opt for the more traditional definition 
9iiv,s = e gitv,E, which implies no factor of g s in the prefactor of W. 
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2.2.2 a' corrections 

The Kahler potential receives corrections at each order in the a' expansion. In the effec- 
tive supergravity description these correspond to higher derivative terms. The leading a' 
contribution comes from the i? 4 term and it leads to the following Kahler potential [ppf] : 

h 2 In | V + -4r7K I ~ -21nV - -JL-. (2.9) 



Here £ is given by £ = — 2*2^ > where x = 2 (/ii,i — ^2,1) is the CY Euler number, and the 

Riemann zeta function is £(3) ~ 1.2. Denoting for convenience £ = £/<7s^ 2 , eq. (|2.S| ) implies 
to leading order the following contribution to V (again, up to a prefactor containing powers 
of g s and Mp): 

(| 2 + 7£V + V 2 ) otw2 

*v> = ^7^7 — ^ ^ zife' (2 - 10) 

'V-£) (2V + £) 4 S« V 3 



where V 3> £ 3> 1 in order for perturbation theory to be valid. 
2.2.3 g s corrections 

The Kahler potential receives also string loop corrections. At present, there is no explicit 
derivation of these corrections from string scattering amplitudes for a generic CY compact- 
ification. Nevertheless, it has been possible to conjecture their form indirectly [20]: 



where 



and 



hi 



^g s Cf K (U,U)(a u t l ) , s 

SK 5* ~ E y > ( 2 - 12 ) 

1=1 



In (2.11), 5K^^ comes from the exchange of closed strings, carrying Kaluza-Klein momen- 



tum, between D7- and D3-branes. The expression ( |2.12j ) is valid for vanishing open string 
scalars and is based on the assumption that all the h\ \ 4-cycles of the CY are wrapped by 
D7-branes. The other term SKW \ in ( |2 . 1 2| ) is due, from the closed string perspective, to 
the exchange of winding strings between intersecting stacks of D7-branes. 

In addition, in ( |2.12| ) the linear combination [ant 1 ) of the volumes of the basis 2- 
cycles is transverse to the 4-cycle wrapped by the i-ih D7-brane, whereas in (|2.13j ) it gives 
the 2-cycle where two D7-branes intersect. The functions C^ K (U, U) and C^{U, U) are, in 
principle, unknown. However, for our purposes they can be viewed as 0(1) constants 4 since 



4 In the T 6 / (Z2 x Z2) orientifold case, where these constants can be computed explicitly |21| , they turn 
out to be, in our conventions, of 0(1) for natural values of the complex structure moduli: Re((7) ~ lm(U) ~ 
0{1). Note that Q uses conventions different from ours. 
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the complex structure moduli are already stabilised at the classical level by background 
fluxes. 

Comparing ( 2.12| ) with (|2.9| ), we notice that 5K^^ is generically leading with respect 



to 5K( a r^ due to the presence of the linear combination (ant ) in the numerator of ( 2.12] ) 



and the fact that each 4-cycle volume Tj = g&ijfc^i has to be fixed larger than the string 
scale in order to trust the effective field theory. However, in ref. p2j it was discovered 
that for an arbitrary CY background, the leading contribution of ( 2.12| ) to the scalar 
potential is vanishing, so leading to an extended no-scale structure. This result renders 
8Vr g A generically subleading with respect to <5V( a /) and the final formula can be expressed 
as (neglecting a prefactors with powers of g s and Mp): 



5V, 



l—loop 



(gs) 



{g s C« K ) 2 a ikaij K° k3 -2SK 



w 

(gs) 



1 W„ 2 



V 2 ' 



(2.14) 



Notice that for branes wrapped only around the basis 4-cycles the combination appearing in 
the first term degenerates to a^a^K^,- = K~. The fact that $V( gs ) is generically subleading 
with respect to 5V( a /) can be easily seen by recalling the generic expression of the tree-level 
Kahler metric for an arbitrary CY: K% = t i t- j /(2V 2 ) - (k i ] k t k )~ 1 /V ~ l/(Vt). Therefore 
the ratio between 5Vr a i\ given by ( |2.10D (with £ ~ C(l) for known CY three-folds) and the 



expression ( [2.14 ) for 8Vr g \, scales as $Vr a i)/5Vr gs -\ ~ g s 3 ^ 2 t 3> 1 due to the fact that the 
size of each two-cycle has to be fixed at a value larger than 1 (in string units) in order to 
trust the effective field theory and, in addition, the string coupling has to be smaller than 
unity in order to be in the perturbative regime. 



2.3 Moduli stabilisation in LVS 



Combining (|2.§|), ( |2.10 ) and ( p. 14 ), we obtain the final form of the Kahler moduli effective 
scalar potential, that contains all the leading order corrections to the vanishing tree-level 
part (with all prefactors included): 

V = Vtree + 6V [np) + + SV^f* = 



g s e K ^M P 



+ 



8vrV 2 

3/2 

,9s 



K^ajAj a i A i e-( a ^ +aiTi ) + 4W ^ a^n cos(a i 6 i ) e - as 



+ £(«(^) 2 (^-a*)-^ 



i 

w 



4V 



(2.15) 



where Aij = ^ = kij k t k . Ref. |^] derived the topological conditions that an arbitrary CY 
has to satisfy in order for the general potential ( 2.15| ) to have a non-supersymmetric AdS 
minimum at exponentially large volume. These conditions can be summarised as follows: 

1. h 2) i > h 1} i > 1 <^ £ ~ {h 2 ,i - /ii,i) > 0, 

2. The CY 3- fold has to have at least one Kahler modulus corresponding to a blow-up 
mode resolving a point-like singularity (the volume of a del Pezzo 4-cycle). 
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These conditions give rise to two different LVS. In the first case, the CY has a typical 
Swiss-cheese topological structure where there is just one LARGE cycle controlling the 
size of the overall volume and all the other 4-cycles are blow-up modes. The string loop 
corrections turn out to be negligible. On the other hand, in the case of fibred CY manifolds, 
with the presence of modes which do not resolve point-like singularities or correspond to 
the overall volume modulus, g s corrections play a crucial role to lift the fibration moduli 
that are left unfixed by 8Vr np \ +5V( a >} (more precisely, the string loop corrections for these 
modes are always dominant compared to non-perturbative effects). Let us now describe 
these two cases in more detail. 

2.3.1 Swiss-cheese Calabi-Yaus 

The original example of LVS of [|| , and the one we will study in most detail, is the degree 
18 hypersurface in CP^ 1 1 6 9] wnose volume is given by 

v = ^( T ' /2 - T - /2 )- < 2 - 16 > 

where tj, and t s are the two Kahler moduli and the subscripts b and s stand for big and 



small respectively. The general expression (2.15) for the scalar potential, in this case takes 
the form 

g„e K «M* /\ e - 2a °T° T s e~ a ^ v a p \ 

v = ^^[ x ^^^-»^^ + V3 + ^ s + v^)> (2 - 17) 

with 

X = 2AV2a 2 s Al p = 4a s A s W , v = ^|f, a = (g s Cf K W ) 2 , p= (gfgglg! , 

±gJ (9v/2) 1 

(2.18) 

Also, in ( 2.17] ) the minimisation with respect to the axion b s has already been performed. 



For natural values of the tree-level superpotential Wq ~ 0(1), the scalar potential ( 2.17| ) 



admits a non-supersymmetric AdS minimum at exponentially large volume due to the 
interplay of a' and non-perturbative effects. This minimum is located at 

V ~ W e asTa > t s ~ i 2/3 > 1 . (2.19) 

The string loop corrections can be safely neglected since they are subdominant relative to 
the other corrections due to inverse powers of V and factors of g s . There are several ways 
to up-lift this minimum to Minkowski or dS: adding D3 branes ||, considering D-terms 
from magnetised D7 branes [27] or F-terms from a hidden sector [^] etc. 

An immediate generalisation of the CP^ 1 1 6 9] m °del is given by the so called 'Swiss- 
cheese' Calabi-Yaus, whose volume looks like 

small \ 

T b /2 ~ E X ^'j ' a > °' ^ > V* = 1, ...,N small . (2.20) 
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Examples having this form with h\ : \ = 3 are the the degree 15 hypersurface embedded in 
CP^ 3335] and the degree 30 hypersurface in CPA x 3 10 15] \—\' More generally, in ref. [25] 
it was proved that examples of Swiss-cheese CY 3-folds with h\ t i = n + 2, 0<n<8, can 
be obtained by starting from elliptically fibred CY manifolds over a del Pezzo dP n base 5 , 
and then performing particular flop transitions that flop away all n CP 1 -cycles in the base. 
In this case, assuming that all the small cycles get non-perturbative effects, the 4-cycle t&, 

3/2 

controlling the overall size of the CY, is stabilised exponentially large, V — ar b ~ Woe aiTi , 
while the various 4-cycles, Tj, controlling the size of the 'holes' of the Swiss-cheese, get fixed 
at small values Tj ~ O(10), Vi = 1, N sma u. However, in ref. [24] it was discovered that 
the Swiss-cheese structure of the volume is not enough to guarantee that all the rigid 'small' 
cycles Tj can indeed be stabilised small. In fact, a further condition is that each rigid 'small' 



cycle Tj must be del Pezzo. In [24|, there are 3 examples of Swiss-cheese CY 3-folds with 
h\ t \ = 4 where just one 4-cycle has the topology CP 2 (and so it is dPo). 

2.3.2 Fibred Calabi-Yaus 

The first examples of LVS with a topological structure more complicated than the Swiss- 
cheese one, were discovered in Q. The authors focused on a K3 fibred CY with /i2,i > 
h± : i = 3, obtained by adding a blow-up mode to the geometry CP^ 1 2 2 6 ]. The volume 
reads: 



V 



a 



0TT2 - 7T, 



3/2 



iiTx - cryr. 



3/2 



(2.21) 



where the constants a and 7 are positive, and ii is the volume of the CP 1 base of the K3 
fibration. Working in the parameter regime T2 > T\ 3> t s , where the volume of the CY 
is large, while the blow-up cycle t s remains comparatively small 6 , the general expression 
( 2.15 ) for V becomes (having already minimised V with respect to the axion b s = ImT s ): 

2a s r 3 Tse -**T S v f A B Ct \ W 2 



V 



8vr 



V 



V 2 



+ 7« + -o 



with A, /x and v as given in ( 2.18j ) and 
A 



(9sCi 



KK\2 



D 



V 3 



iaC w 



12 ! 



C = 2(ag s C, 



KK\2 
2 ) ■ 



(2.22) 



(2.23) 



9\/2a7 ' 

It is evident that the leading <5V( a n + 8Vt np \ part of the potential depends only on two 
Kahler moduli, V and t s , instead of all three. And, in fact, it turns out to be of exactly the 
same form as ( p. 17 ) above. Hence, viewing V, r s and t\ as the three independent moduli 
(instead of t±, T2 and r s ), it is clear that, without taking into account the subleading g s 
corrections, n is a flat direction of the scalar potential. Also, it is evident that at this 



order, V and t s are stabilised as in the CP^ 1169 j model of Subsection 2.3.1 

2/3 

(rs) 



2 07 



and 



(V) 



3 ocy 



Aa s A s 



,a s (T s ) 



(2.24) 



5 A del Pezzo dP n surface is obtained by blowing-up CP 2 (or CP 1 x CP 1 ) on < n < 8 points. 
6 in this limit t\ ~ TilsfrL > \Jt\-, corresponding to interesting geometries having the two dimensions of 
the base, spanned by the cycle t\, larger than the other four of the K3 fibre, spanned by n. 
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9s 


e 


Wo 


a3 


^3 


a 


7 




C« K 


°12 


<r.> 


(n) 


(V) 


LV 
SV 


0.1 
0.3 


0.4 
0.9 


1 
100 


tt/5 


1 
1 


0.5 
0.13 


0.39 
3.65 


0.1 
0.15 


0.1 

0.08 


5 
1 


10.5 
4.3 


10 6 
9 


3 • 10 13 
1710 



Table 1: Some model parameters. 



Obviously, loop corrections shift insignificantly the VEVs of these two moduli. However, 
g s corrections are crucial to generate a potential for n that admits a minimum at 



-3/2 \8AVJ 



. r> : , 32AC 
l + ( sl gnS) A /l + -^- 



(2.25) 



Some concrete numerical choices for the various underlying parameters, without any 
fine-tuning, are listed in Table 1. The 'LV case gives very large volumes, V — 10 13 and 
the modulus t\ is stabilised at hierarchically large values, T2 > T\ ^> t s . The string scale 
and the gravitino mass turn out to be 

M. = -^~10"GeV, m 3/2 = = ^""^ ~ IQTeV. (2.26) 

This gives a solution of the hierarchy problem but the huge value of the volume destroys the 
standard picture of gauge coupling unification. The 'SV case instead has V ~ 10 3 much 
smaller (and so with M s ~ Mqut and a very high gravitino mass). This set of parameter 
values is not chosen just in relation to GUT theories but also in order to provide observable 
density fluctuations for the inflationary model of jR]]. In that model, the inflaton is the 
modulus T\ , whose potential is loop-generated, and the main feature of the model is that 
it produces detectable gravity waves. 



More general examples of this kind of LVS have been discovered in 25]. These authors 
noticed that starting from an elliptically fibred CY over a dP n base, and then flopping away 
only r < n (instead of all n) of the CP 1 -cycles in the base, one obtains another elliptically 
fibred CY (instead of a Swiss-cheese one), whose volume looks like: 

r 

V = Vol{X r ^ r )-Y,^'\ \ > Vi = l,...,r, (2.27) 

i=l 

where V„_ r is the resulting elliptical fibration over a dP n - r base. It is natural to expect 
that the scalar potential for these examples has an AdS minimum at exponentially large 
volume, together with (h\ t i — N sma u — 1) = n — r flat directions that will be lifted by g s 
corrections. 

We should note that string loop corrections can play an important role for compacti- 
fications on Swiss-cheese CY manifolds as well. Namely, they can be crucial, even in this 
case, to achieve full moduli stabilisation when the topological condition, that all rigid 4- 
cycles be del Pezzo, is not satisfied or when one imposes the phenomenological condition 
that the 4-cycles supporting chiral matter do not get non-perturbative effects [ 23 1 . ' 



7 Also D-terms could play a significant role as pointed out still in 
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3. Effective potential at finite temperature 



At nonzero temperature, the effective potential receives a temperature-dependent contri- 
bution. The latter is determined by the particle species that are in thermal equilibrium 
and, more precisely, by their masses and couplings. In this Section, we review the general 
form of the finite temperature effective potential and discuss in detail the establishment of 
thermal equilibrium in an expanding Universe. In particular, we elaborate on the relevant 
interactions at the microscopic level. This lays the foundation for the explicit computation, 



in Section |6.l| , of the finite temperature effective potential in LVS. 
3.1 General form of temperature corrections 

The general structure of the effective scalar potential is the following one: 

Vtot = V + Vt, (3.1) 

where Vo is the T = potential and Vt the thermal correction. As discussed in Section 0, 
Vo has the general form: 

V = 5V {np) +6V {a , ) + 5V {gs) , (3.2) 

where the tree level part is null due to the no-scale structure (recall that we are studying 
the scalar potential for the Kahler moduli), 5V( np ) arises due to non-perturbative effects, 
<5V( Q /) are a' corrections and the contribution 5V( 9s ) comes from string loops and, as noticed 



in [22], matches the Coleman- Weinberg potential of the effective field theory. In addition, 



(5V( 9s ) has an extended no-scale structure, which is crucial for the robustness of LVS since 
it renders 5V( 9s ) subleading with respect to 5V( np ) and 5V( a >y 

On the other hand, the finite temperature corrections Vt have the generic loop expan- 
sion: 

V T = V^' loop + v*- loops + ... . (3.3) 

The first term V^, loop is a 1-loop thermal correction describing an ideal gas of non- 
interacting particles. It has been derived for a renormalisable field theory in flat space 
in [33], using the zero-temperature functional integral method of [34], and reads: 



yl-loOp = ± 



T 
2^2 



dx x 2 In (l T e -V^)W) , (3.4) 



where the upper (lower) signs are for bosons (fermions) and m is the background field 
dependent mass parameter. At temperatures much higher than the mass of the particles 



in the thermal bath, T 3> m(ip), the 1-loop finite temperature correction (3^) has the 
following expansion: 

^=-?« + ^ir +0 Wrf). (3-5) 

where for bosons a = 1 and for fermions a = 7/8. The generalisation of (3.5) to super- 
gravity, coupled to an arbitrary number of chiral superfields, takes the form [31]: 

V T ' loop = ~^ (g B + Iqf) + Ya ( TrM * + TrM /) + ° (™') ' ( 3 ' 6 ) 
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where gp and gp are, respectively, the numbers of bosonic and fermionic degrees of freedom 
and Mf, and Mf are the moduli-dependent bosonic and fermionic mass matrices of all the 
particles forming the thermal plasma. 

If the particles in the thermal bath interact among themselves, we need to go beyond 
the ideal gas approximation. The effect of the interactions is taken into account by eval- 
uating higher thermal loops. The high temperature expansion of the 2-loop contribution 
looks like: 

V^ loops = a 2 T* fi(9i)\ + P2T 2 (TrMl + TrM 2 ) ^£ ffa)} + ... , (3.7) 

where «2 and are known constants, i runs over all the interactions through which 
different species reach thermal equilibrium, and the functions fi are determined by the 
couplings gi and the number of bosonic and fermionic degrees of freedom. For example, 
for gauge interactions f(g) = const x g 2 , whereas for the scalar A</> 4 theory one has that 
/(A) = const x A. 

Now, since we are interested in the moduli-dependence of the finite temperature correc- 
tions to the scalar potential, we can drop the first term on the RHS of ( |3.6[) and focus only 
on the T 2 term, which indeed inherits moduli-dependence from the bosonic and fermionic 
mass matrices. However, notice that in string theory the various couplings are generically 



functions of the moduli. Thus, also the first term on the RHS of (3.7) depends on the 
moduli and, even though it is a 2-loop effect, it could compete with the second term on the 
RHS of Q3.6| ), because it scales as T 4 whereas the latter one scales only as T 2 . This issue 
has to be addressed on a case by case basis, by studying carefully what particles form the 
thermal bath. 

3.2 Thermal equilibrium 

In an expanding Universe, a particle species is in equilibrium with the thermal bath if 
its interaction rate, T, with the particles in that bath is larger than the expansion rate 
of the Universe. The latter is given by H ~ g^ 2 T 2 /Mp, during the radiation dominated 
epoch, with g* being the total number of degrees of freedom. Thermal equilibrium can 
be established and maintained by 2 *-* 2 interactions, like scattering or annihilation and 
the inverse pair production processes, and also by 1 *-* 2 processes, like decays and inverse 
decays (single particle productions). Let us now consider each of these two cases in detail. 

3.2.1 2^2 interactions 

In this case the thermally averaged interaction rate can be inferred on dimensional grounds 
by noticing that: 

<r>~^, (3.8) 

where (t c ) is the mean time between two collisions (interactions). Moreover 

t c ~ — , (3.9) 
no~v 
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where n is the number density of the species, a is the effective cross section and v is the 
relative velocity between the particles. Thus (r) ~ n(av). For relativistic particles, one 
has that (v ) ~ c (= 1 in our units) and also n ~ T 3 . Therefore 

<r) ~ (a)T 3 . (3.10) 

The cross-section a has dimension of (length) 2 and for 2 <-> 2 processes its thermal average 
scales with the temperature as: 

1. For renormalisable interactions: 

T 2 

(a) ~ a 2 (3.11) 

X ' (T 2 + M 2 ) 2 V ' 

where a = <7 2 /(47r) (g is the gauge coupling) and M is the mass of the particle 
mediating the interactions under consideration. 

a) For long-range interactions M = and ( p. 11 ) reduces to: 

(a) ~ a 2 T" 2 (r) ~ a 2 T. (3.12) 

This is also the form that (3.11) takes for short-range interactions at energies 
E » M. 

b) For short-range interactions at scales lower than the mass of the mediator, the 
coupling constant becomes dimensionful and ( 3.11| ) looks like: 

(^)~« 2 ^ * {T) ~ a2 w*- (3 - 13) 

2. For processes including gravity: 

a) Processes with two gravitational vertices: 

(a) ~ (r) ~ d^x, (3.14) 

Mp w Mp' v 7 

where d is a dimensionless moduli-dependent constant. 

b) Processes with one renormalisable and one gravitational vertex: 



where (i is the same moduli-dependent constant as before. 

Let us now compare these interaction rates with the expansion rate of the Universe, 
H ~ gl^ 2 T 2 /Mp, in order to determine at what temperatures various particle species reach 
or drop out of thermal equilibrium, depending on the degree of efficiency of the relevant 
interactions. 
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Figure 2: QCD scattering process qq — > gg through which quarks and gluons reach thermal 
equilibrium. 



l.a) Renormalisable interactions with massless mediators: 

(r) > H a 2 T > gl /2 T 2 M p l T < a 2 g~ 1/2 M P . (3.16) 

QCD processes, like the ones shown in Figure |2], are the main examples of this kind 
of interactions. The same behaviour of a is expected also for the other MSSM gauge 
groups for energies above the EW symmetry breaking scale. Therefore, MSSM parti- 
cles form a thermal bath via strong interactions for temperatures T < a^g* Mp ~ 



10 15 GeV [36]. 



l.b) Renormalisable interactions with massive mediators: 

lT)> H « ^>£'% - <T<M. ,3.7, 

Examples of interactions with effective dimensionful couplings are weak interactions 
below Mew- I n this case, the theory is well described by the Fermi Lagrangian. An 
interaction between electrons and neutrinos, like the one shown in Figure ||, gives rise 
to a cross-section of the form of ( |3.13| ) : 

~ ~ ik T2 ' (3 ' I8) 

where a w is the weak fine structure constant and p ~ T. Thus, neutrinos are coupled 
to the thermal bath if and only if 

/ 1/2 , f4 \ V3 

T > 9 Ayi^ ~ 1 MeV. (3.19) 



2. Gravitational interactions: 

a) { T)>H « ^> 3 l' 2 ^ - T> g l"'^. (3.20) 
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Figure 3: Weak interaction between electrons and neutrinos through which they reach thermal 



equilibrium. 




Figure 4: Scattering process $g — ► gg through which the modulus $ and gluons can reach thermal 
equilibrium. 

i \ ,r\^tr ^ f7i3 2T '' 1/2 T 2 gl /2 M P 

b) (T)>H o Vd 1 ^> 9 * — =* T>-^ T . (3.21) 

As before, case (a) refers to 2 «-> 2 processes with two gravitational vertices, whereas 
in case (b) one vertex is gravitational and the other one is a renormalisable interac- 
tion. A typical Kahler modulus of string compactifications generically couples to the 
gauge bosons of the field theory, that lives on the stack of branes wrapping the cycle 
whose volume is given by that modulus. Scattering processes, annihilation and pair 
production reactions, that arise due to that coupling, all have cross-sections of the 
form ( |3. 14 ) and ( 3.15| ). For all the Kahler moduli in KKLT constructions d ~ O(l) 



and so (r) is never greater than H for temperatures below the Planck scale, for both 
cases (a) and (b). Therefore, those moduli will never thermalise through 2 «-> 2 pro- 
cesses. However, we shall see in Section || that the situation is different for the small 
modulus in LVS, since in that case d ~ V 2 S> 1. A typical 2 <-> 2 process of type 
(b), with a modulus $ and a non-abelian gauge boson g going to two g's, is shown 
in Figure ||. Here $ denotes the canonically normalized field, which at leading order 
in the large-volume expansion corresponds to the small modulus. We will give the 
precise definition of <3? in Section ||. 



3.2.2 1 <-> 2 interactions 

In order to work out the temperature dependence of the interaction rate for decay and 
inverse decay processes, recall that the rest frame decay rate 1^ does not depend on 
the temperature. For renormalisable interactions with massless mediators or mediated by 
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particles with mass M at temperatures T > M, it takes the form: 

~ am, (3.22) 

where m is the mass of the decaying particle and a ~ g 2 , with g being either a gauge or a 
Yukawa coupling. On the other hand, for gravitational interactions or for renormalisable 
interactions mediated by particles with mass M at temperatures T < M, we have (M = Mp 
in the case of gravity): 

rg> ~ D £, (3.23) 

with D a dimensionless constant (note that in the case of gravity D = yd, where d is the 
same moduli-dependent constant as in Subsection 3.2.1| ). 



( Ft) 

Now, the decay rate that has to be compared with H is not F D , but its thermal 
average (Tp). In order to evaluate this quantity, we need to switch to the 'laboratory 
frame' where: 

i m 
~E 

and then take the thermal average: 



T D =T^V^ = T^, (3.24) 



(T D )=T^^. (3.25) 

In the relativistic regime, T > m, the Lorentz factor 7 = (E)/m ~ T/m, whereas in the 
non-relativistic regime, T < m, 7 = (E)/m ~ 1. 

Notice that, by definition, in a thermal bath the decay rate of the direct process is 
equal to the decay rate of the inverse process. However, for T < m the energy of the final 
states of the decay process is of order T, which means that the final states do not have 
enough energy to re-create the decaying particle. So the rate for the inverse decay, Tip, is 
Boltzmann-suppressed: Tjp ~ e~ m l T . Hence, the conclusion is that, for T < m, one can 
never have Tp = Tfp and thermal equilibrium will not be attained. Let us now summarize 
the various decay and inverse decay rates: 

1. Renormalisable interactions with massless mediators or mediated by particles with 
mass M at T > M: 

<r D > = (3.20) 

g m, tor 1 ^ m , 



g 2 ^, for T>m 
g 2 m (f) 3/2 e~ m / T , for T < w . 



<pid)*< a Jur:;; r (3.27) 



Therefore, particles will reach thermal equilibrium via decay and inverse decay pro- 
cesses if and only if 



g* 
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2. Gravity or renormalisable interactions mediated by particles with mass M at T < M: 

(T D )^\ D ^ fOTT ^ (3.29) 

(/D> \D^m V2 e- m/T , forT< m ^ 

with M = Mp in the case of gravity. Therefore, particles will reach thermal equilib- 
rium via decay and inverse decay processes if and only if 



In the case of gravitational interactions, ( 3.31] ) becomes 

1/3 




(3.32) 



In KKLT constructions, D ~ 0(1) and m ~ TO3/2. So ( |3.32 ) can never be satisfied 



and hence moduli cannot reach thermal equilibrium via decay and inverse decay 
processes. However, we shall see in Section [5] that in LVS one has D~ V> 1 and so 
1 *-* 2 processes could, in principle, play a role in maintaining thermal equilibrium 
between moduli and ordinary MSSM particles. 

4. Moduli masses and couplings 

As we have seen in the previous section, the temperature, at which a thermal bath is 
established or some particles drop out of thermal equilibrium, depends on the masses and 
couplings of the particles. To determine the latter, one needs to use canonically normalised 
fields. In this Section, we study the canonical normalisation of the Kahler moduli kinetic 
terms and use the results to compute the masses of those moduli and their couplings to 
visible sector particles. 

4.1 Single-hole Swiss-cheese 

We start by focusing on the simplest Calabi-Yau realisation of LVS, the 'single-hole Swiss- 
cheese' case described in Subsection |2.3.1| (i.e., the degree 18 hypersurface embedded in 

CP [i.i,i,6,9])- First of a11 ' 

we shall review the canonical normalisation derived in |]l7| . In 
order to obtain the Lagrangian in the vicinity of the zero temperature vacuum, one expands 
the moduli fields around the T = minimum: 

n = (n) + 5n , (4.1) 
t s = (t s ) + 5t s . (4.2) 
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where (t&) and (r s ) denote the VEV of n and t s . One then finds: 

£ = K^Sr^i&ri) - (Vb) - ^t^t,- + 0(<5r 3 ) , (4.3) 

where i = b,s and (Vq) denotes the value of the zero temperature potential at the minimum. 
To find the canonically normalized fields $ and x> let us write Sn and St s as: 

Then the conditions for the Lagrangian ( [4.3| ) to take the canonical form: 

C = \d^d^ + \d^X - {Vo) - \m%$> 2 - X -m 2 xX 2 (4.5) 



are the following: 

K ij(va)i{vp)j = 6 a p and ^yCj{v a )i{yp)j = m 2 a 5 af} . ( 



These relations are satisfied when v$, v x (properly normalised according to the first of 
(|4.6| )) and m|, mi are, respectively, the eigenvectors and the eigenvalues of the mass- 
squared matrix {M 2 ) tj = \ {K' 1 ) fk V- kj . 

Substituting the results of [17| for v§ and v x in (fO|), we can write the original Kahler 
moduli Sri as (for a s T s S> 1): 



Sn 



V6(n)^(T s )^) * + ( J\(n) )j,-o (vV") $ + o (v 2 / 3 ) x , (4.7) 



x/2 



Sr. = ( 2 -f infers) 



1/4 )^ + (^)^^( Vl/2 ) $ + 0(1) - (4 ' 



As expected, these relations show that there is a mixing of the original fields. Nevertheless, 
Sn is mostly x an d St s is mostly $. On the other hand, the mass-squareds are [17|: 

2 m /,, 2 n fg s e K ™\ 24V2i/ a 2(r s ) 1 /2 /lnV\ 2 2 ,„ , 

m| ~ Tr (M ) ~ ^- g/ — M| ~ — M| (4.9) 



8vr J V 2 V V , 

2 Det(M 2 ) {g s e K <*\ 27u 2 Mp . , 



x Tr (M 2 ) \ 8vr J 4a s (r s )V 3 ^ V 3 lnV' 

We can see that there is a large hierarchy of masses among the two particles, with $ being 
heavier than the gravitino mass (recall that m 3 / 2 ~ Mp/V) and \ lighter by a factor of \fV . 



Using the above results and assuming that the MSSM is built via magnetised D7 branes 
wrapped around the small cycle, we can compute the couplings of the Kahler moduli fields 
of the CP^ model to visible gauge and matter fields. This is achieved by expanding 

the kinetic and mass terms of the MSSM particles around the moduli VEVs. The details 
are provided in Appendix [A], where we focus on T > Mew since we are interested in 
thermal corrections at high temperatures. This, in particular, means that all fermions and 
gauge bosons are massless and the mixing of the Higgsinos with the EW gauginos, that 
gives neutralinos and charginos, is not present. We summarise the results for the moduli 
couplings in Tables |2| and ||. 
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Gauge bosons (F^F^) 


Gauginos (AA) 


Matter fermions (ipip) 


Higgsinos (HH) 


X 


1 


1 


No coupling 


l 


M p InV 


VlnV 


VlnV 




Vv 


1 


No coupling 


1 




V 3 / 2 lnV 


v^lnV 



Table 2: CP^ i i 6 9] case: moduli couplings to spin 1 and 1/2 MSSM particles for T > Mew- 





Higgs (HH) 


Higgs-Fermions (Hipip) 


SUSY scalars (ip<p) 


x 2 


$2 


X 


M P 


1 


M P 


Mp 


M P 


V 2 (lnV) 2 


MpVV 3 


V 2 (lnV) 2 


~w 


~W 




M P 


1 


M P 


Mp 


M P 


V s / 2 (InV) 2 


M P V 5 /6 


V 5 /2(l n v) 2 


y&72 





Table 3: CPS 1 16 9 ] case: moduli couplings to spin and 1/2 MSSM particles and cubic self- 
couplings for T > Mew- 



4.2 Multiple-hole Swiss-cheese 

Let us now consider the more general Swiss-cheese CY three-folds with more than one 



small modulus and with volume given by (2.20). In this case we find: 



9 
4 



(4.11) 



— <iC 1 and also we have kept only the leading (in the limit r b >5> 73 Vi) 



where q = 

contribution in each term. Notice that the mixed terms are subleading compared to the 



diagonal ones. So, to start with, one can keep only the first line in ( 4.11 ). Then at leading 
order the canonically normalized fields x an d i = 1, ...,iV sma (/, are defined via: 



5r b 



5n 



/3l 



■.(n^in^i-ofv 1 / 2 )^. (4.12) 



As was to be expected, this scaling with the volume agrees with the behaviour of 5r b and 
5t s in ( f4.7|) , (4.8). Now, let us work out the volume scaling of the subdominant mixing 
terms since it is important for the computation of the various moduli couplings. Proceeding 
order by order in a large-V expansion, we end up with: 



5T b ~o(v 2 / 3 )x + Y,°( vl/6 )® 



8n ~ O (V 1/2 j $i + O (1) x + Yl ( v ~ 1/2 ) $3 



(4.13) 
(4.14) 



This shows that the mixing between the small moduli is strongly suppressed by inverse 
powers of the overall volume, in accord with the subleading behaviour of the last term 
in ( f4,ll| ). Furthermore, the fact that the leading order volume-scaling of ( |4.13| )-( ^.14 ) 
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is the same as fl4.7|) -( f4T8|) , implies that all small moduli behave in the same way as the 
only small modulus of the CP^ 1 1 6 g , model. Hence, if all the small moduli are stabilised 
by non-perturbative effects, the moduli mass spectrum in the general case will look like 
fl4,9| )-( |4~Tc| ), with Q4.9| ) valid for all the small moduli. In addition, if we assume that all 
the 4-cycles corresponding to small moduli are wrapped by MSSM D7 branes, the moduli 
couplings to matter fields are again given by Tables ^ and || where now <3? stands for any 
small modulus 

However, in general the situation may be more complicated. In fact, the authors of 
|23|1 pointed out that 4-cycles supporting MSSM chiral matter cannot always get non- 
perturbative effects. 8 A possible way to stabilise these 4-cycles is to use g s corrections 
as proposed in B. In this case, the leading-order behaviour of Q4.9|) should not change: 



m % ~ Tjf - 9 However, the moduli couplings to MSSM particles depend on the underlying 
brane set-up. So let us consider the following main cases: 

1. All the small 4-cycles are wrapped by MSSM D7 branes except r np which is respon- 
sible for non-perturbative effects, being wrapped by an ED3 brane. It follows that 
the MSSM couplings of <& np are significantly suppressed compared to the MSSM cou- 
plings of the other small cycles (still given by Tables [2] and |3|). This is due to the 
mixing term in ( [4.14 ) being highly suppressed by inverse powers of V. 



2. All the small 4-cycles are wrapped by MSSM D7 branes except r np which is supporting 
a pure SU(N) hidden sector that gives rise to gaugino condensation. This implies 
that the coupling of & n p to hidden sector gauge bosons will have the same volume- 
scaling as the coupling of the other small moduli with visible sector gauge bosons. 
However, the coupling of the MSSM 4-cycles with hidden sector gauge bosons will be 
highly suppressed. 

3. All the small 4-cycles Tj support MSSM D7 branes which are also wrapped around the 
4-cycle responsible for non-perturbative effects T np , but they have chiral intersections 
only on the other small cycles. In this case, the coupling of <J? np to MSSM particles 
would be the same as the other <3?j. However, if r np supports an hidden sector that 
undergoes gaugino condensation, the coupling of the MSSM 4-cycles with the gauge 
bosons of this hidden sector would still be highly suppressed. 

4.3 K3 Fibration 



We turn now to the K3 fibration case described in Section 2.3.2. We shall consider first 



the 'LV case, in which the modulus related to the K3 divisor is fixed at a very large value, 



8 This is because an ED3 wrapped on the same cycle will have, in general, chiral intersections with the 
MSSM branes. Thus the instanton prefactor would be dependent on the VEVs of MSSM fields which are 
set to zero for phenomenological reasons. In the case of gaugino condensation, this non-perturbative effect 
would be killed by the arising of chiral matter. 

9 It may be likely that m% . depends on subleading powers of (In V) due to the fact that the loop corrections 
are subdominant with respect to the non-perturbative ones (see 0|), but the main V -2 dependence should 
persist. 
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and then the 'SV case, in which the overall volume is of the order V ~ 10 3 and the K3 
fiber is small. 

In order to compute the moduli mass spectroscopy and couplings, it suffices to canon- 
ically normalise the fields just in the vicinity of the vacuum. The non-canonical kinetic 
terms look like (with e = \J (t s )/(ti)): 

- ^^(<5t 2 )^(«5t s ) + -J^-d^d^Sn) + -JVL-d^T.WVT.). (4.15) 
2(r 2 )^ r 2(T 2 ) i! ^ 8(r 2 )(r s ) 



Large K3 fiber 

In the 'LV case where the K3 fiber is stabilised at large value, eCl. Therefore at leading 
order in a large volume expansion, where (t 2 ) > (ti) 3> (r s ), all the cross-terms in ( 4.15| ) 
are subdominant to the diagonal ones, and so can be neglected: 

Ckin * 47^2^(^)^(*n) + ^^(5r 2 )^(5r 2 ) + J^— dli (6T,)d>>(6T,). (4.16) 

Therefore, at leading order the canonical normalisation close to the minimum becomes 
rather easy and reads: 

Sn = v / 2(n)xi ~ o (v 2 / 3 ) xi, (4.17) 

<5r 2 = (r 2 ) X2 ~ O (V 2/3 ) X 2, (4.18) 
8r^^^^^0(yy^. (4.19) 

However, in order to derive all the moduli couplings, we need also to work out the leading 
order volume-scaling of the subdominant mixing terms in JDg) and (14191) . This can be 
done order by order in a large-V expansion and, after some algebra, we obtain: 

5r 1 =a 1 {T 1 ) Xl + a 2 y^(T s ) 3 / 2 X2 + a 3 { -^^(T s f^, (4.20) 
( r 2/ V(r 2 ) 

5r 2 = a 4 ^M( Ts) 3/2 Xi + as(T2)x2 + a6 ^/^ (Ts) 3/4 $; (4 21) 

<5r s = a 7 ^(r s ) Xl + a 8 (r s ) X2 + agin) 1 / 4 ^^) 1 / 4 ^ (4.22) 

\ T 2) 

where the Qj, i = 1,...,9 are 0(1) coefficients. The volume-scalings of ( 4.20| ), ( [1.21 ) and 
(4.22) are the following: 



Sn ~ o (v 2/3 J xi + o (y- 1,z ) X 2 + o (v 1/6 J (4.23) 

5r 2 ~ (V 1/3 ) xi + (V 2/3 ) X 2 + O (V 1/6 ) (4.24) 
<5r s ~ O (1) xi + O (1) X2 + O (V 1 ' 2 ) (4.25) 
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This shows that, if we identify each of t\ and ti with the large modulus r& in the Swiss- 
cheese case, (|423|) and (|3|) have the same volume scaling as (Hj), as one might have 



expected. Moreover, the similarity of ( |4,25| ) and ( |4.8| ) shows that also the small moduli in 
the two cases behave in the same way. Therefore, we can conclude that (4.9) is valid also 
for the K3 Fibration case under consideration: 

m$ ~ ( — J M P . (4.26) 

On the other hand, we need to be more careful in the study of the mass spectrum of the 
large moduli t\ and T2. We can work out this 'fine structure', at leading order in a large- V 
expansion, first integrating out t s and then computing the eigenvalues of the matrix. The 
latter are obtained by multiplying the inverse Kahler metric by the Hessian of the potential 
both evaluated at the minimum. The leading order behaviour of the determinant of this 
matrix is: 

Det (K- l d 2 V) ~ T2 ^f^ , with V~0Tr 2 . (4.27) 
Because m£ 2 3> m 2 , 1 , we have at leading order at large volume: 



m 2 X2 ~ Tr {K- l d 2 V) ~ (4.28) 



mt oi - 1 loTW ~ ^ m p ~ "^4- ( 4 -29) 



Identifying t% with T2, 29| ) simplifies to ~ V -10 / 3 , confirming the qualitative expec- 
tation that the t\ direction is systematically lighter than V in the large- V limit. 

Using the results of this Section and assuming that the MSSM branes are wrapped 
around the small cycle 10 , it is easy to repeat the computations of Appendix |A| for the K3 



fibration. Due to the fact that the leading order V-scaling of ( 4.23 )-( 4.25 ) matches that 
of the single-hole Swiss-cheese model, we again find the same couplings as those given in 
Tables ^ and ||, where now \ stands for any of \\ an( i X2- 

Small K3 fiber 

In the 'SV case where the K3 fiber is stabilised at small value, e ~ 1. Therefore at leading 
order in a large volume expansion, where (7$) S> {t\) > (r s ), the first term in ( 4,15| ) 



is dominating the whole kinetic Lagrangian. Hence we conclude that, at leading order, 
the canonical normalisation of 5t\ close to the T = minimum is again given by ( 4.17] ). 
However, now its volume scaling reads: 

8t\ ~ O (1) xi + (subleading mixing terms) . (4.30) 

To proceed order by order in a large volume expansion, note that the third and the sixth 



term in (4.15) are suppressed by just one power of (T2), whereas the second, fourth and fifth 



10 We also ignore the incompatibility between localising non-perturbative effects and the MSSM on the 
same 4-cycle. 
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term are suppressed by two powers of the large modulus. Thus, we obtain the following 
leading order behaviour for the canonical normalisation of the two remaining moduli: 



5r 2 ~ O (V) xi + O (V) X2 + O (V) $, (4.31) 
5t s ~ O (v 1/2 ) Xi + O (v 1/2 ) $ + subleading mixing terms. (4.32) 



Notice that the canonically normalised field xi corresponds to the K3 divisor r\, whereas 
$ is a mixing of n and the blow-up mode t s . Finally X2 is a combination of all the three 
states, and so plays the role of the 'large' field. The moduli mass spectrum will still be 
given by (4.26), ( |4.28| ) and ( |4.2S| ). However now the volume scaling of (4.29) simplifies to 



m 2 ~ V 2 , confirming the qualitative expectation that xx is also a small field with a mass 
of the same order of magnitude of m$. 

The computation of the moduli couplings depends on the localisation of the MSSM 
within the compact CY. As we have seen in Subsection p. 3. 2 , the scalar potential receives 



non-perturbative corrections in the blow-up mode r s . Therefore, in order for the non- 
perturbative contributions to be non- vanishing, the MSSM branes have to wrap either the 
small K3 fiber t\ or the 4-cycle given by the formal sum r s + t\ with chiral intersections 
on t\. In both cases, we cannot immediately read off the moduli couplings from the results 
of Appendix |A]. This is due to the difference of the leading order volume scaling of the 
canonical normalisation between the 'SV case for the K3 fibration and the Swiss-cheese 
scenario. 11 

However, as we shall see in the next Section, in the Swiss-cheese case, the relevant 
interactions through which the small moduli can thermalise, are with the gauge bosons. 



As we shall see in Section 5.3, these interactions will also be the ones that are crucial for 



moduli thermalisation in the K3 fibration case. Therefore, here we shall focus on them 



only. Following the calculations in Subsection A..1 of Appendix |A|, we infer that if only t\ 
is wrapped by MSSM branes, then the coupling of xx with MSSM gauge bosons is of the 
order g ~ 1/Mp without any factor of the overall volume, while the coupling of $ with 
gauge bosons will be more suppressed by inverse powers of V. On the other hand, if both 
T\ and t s are wrapped by MSSM branes, then the couplings of both small moduli with the 
gauge bosons are similar to the ones in the Swiss-cheese case: g ~ y/V/Mp. Moreover, if 
gaugino condensation is taking place in the pure SU(N) theory supported on t s , then both 
Xx and couple to the hidden sector gauge bosons with strength g ~ Vv/Mp. 

We end this Subsection by commenting on K3 fibrations with more than one blow-up 
mode. In such a case, it is possible to localise the MSSM on one of the small blow-up modes 
and the situation is very similar to the one outlined for the multiple-hole Swiss-cheese. The 
only difference is the presence of the extra modulus related to the K3 fiber, which will couple 
to the MSSM gauge bosons with the same strength as the small modulus supporting the 
MSSM. This is because of the particular form of the canonical normalisation, which, for 



example in the case of two blow-up modes r s \ and r S 2, looks like (|4.30|) and (4.31) together 



11 We stress also that presently there is no knowledge of the Kahler metric for chiral matter localised on 
deformable cycles. 
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with: 



&Tsi ~ O [V 1/2 j Xi + O [V 1/2 J $1 + subleading mixing terms, (4.33) 
5t s2 ~ O (v 1/2 ) Xi + O (v 1/2 ) $2 + subleading mixing terms. (4.34) 



4.4 Modulini 

In this Subsection we shall concentrate on the super symmetric partners of the moduli, 
the modulini. More precisely, we will consider the fermionic components of the chiral 
superfields, whose scalar components are the Kahler moduli. The kinetic Lagrangian for 
these modulini reads: ^ 

£kin = , (4.35) 

where the Kahler metric is the same as the one that appears in the kinetic terms of the 
Kahler moduli. Therefore, the canonical normalisation of the modulini takes exactly the 
same form as the canonical normalisation of the corresponding moduli. For example, in 
the single-hole Swiss-cheese case, we have: 

Sh = (VEWW) ^+(\[%))^~ O * + O (V/ 3 ) x, (4-36) 




Sfs = ( ^#(T5) 3 / 4 (r s )V4 ) ^ + ( ^ ] _| „ o ( ) * + O (1) x. (4.37) 

We focus now on the modulini mass spectrum. We recall that in LVS the minimum is 
non-supersymmetric, and so the Goldstino is eaten by the gravitino via the super-Higgs 
effect. The Goldstino is the supersymmetric partner of the scalar field, which is responsible 
for SUSY breaking. In our case this is the modulus related to the overall volume of the 
Calabi-Yau, as can be checked by studying the order of magnitude of the various F-terms. 
Therefore, the volume modulino is the Goldstino. More precisely, in the CR-j 1 1 6 9 j case, 
X is eaten by the gravitino, whereas the mass of <£> can be derived as follows: 

ml = TrMj = (e G K^V^ + ^)(VjG m + ^p)>, (4.38) 

where the function G = K + ln|VF| 2 is the supergravity Kahler invariant potential, and 
VtGj = Gij — T\jGi, with the connection T\j = K lm diKj m . Equation ( [4.38| ) at leading 
order in a large volume expansion, can be approximated as 



mf ^<e G |(i^(V s G s • ^! (1.39) 



where V S G S ~ G ss — T S SS G S and F s ss ~ K ss d s K sS . In the single-hole Swiss-cheese case, for 
a s r s 3> 1, we obtain: 

™1 - ( - P 36ajA 2 s T s e- 2a ° T ° - s - ^T s e~ a ^ + ^ >. (4.40) 
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Evaluating ( 4.40| ) at the minimum, we find that the mass of the modulino is of the same 



order of magnitude as the mass of its supersymmetric partner <I>: 

m |,W2M^(!f) 2 M^ ra |. (4.4!) 

Similarly, it can be checked that, in the general case of multiple-hole Swiss-cheese Calabi- 
Yaus and K3 fibrations, the masses of the modulini also keep being of the same order of 
magnitude as the masses of the corresponding supersymmetric partners. 

We now turn to the computation of the modulini couplings. In fact, we are interested 
only in the modulino-gaugino-gauge boson coupling since, as we shall see in Section [5], this 
is the relevant interaction through which the modulini reach thermal equilibrium with the 
MSSM thermal bath. This coupling can be worked out by recalling that the small modulus 



t s couples to gauge bosons X as (see appendix |A.l|) : 



{-■gauge ~ , , F^ V F^ . (4.42) 
IVl p 

The supersymmetric completion of this interaction term contains the following modulino- 
gaugino-gauge boson coupling: 

C-^-WF^. (4.43) 
Now, expanding f s around its minimum and going to the canonically normalised fields G 



and A defined as (see appendices A.l and |A7 



Gfj,L> = y/ffiF^ , A = y^A' , (4.44) 

we obtain: 

Hence, by means of fl4.37j ), we end up with the following dimensionful couplings: 

^XX~(t7^)x^XG^, (4.46) 




£**x ~ to^xG^ . (4.47) 



5. Study of moduli thermalisation 



Using the general discussion of Section 3.2 and the explicit expressions for the moduli 



masses and couplings of Section we can now study in detail which particles form the 
thermal bath. Consequently, we will be able to write down the general form that the finite 



temperature corrections of Section 3.1 take in the LVS. 

We shall start by focusing on the simple geometry CPm 1 1 6 g i , and then extend our 
analysis to more general Swiss-cheese and fibred CY manifolds. We will show below that, 
unlike previous expectations in the literature, the moduli corresponding to small cycles 
that support chiral matter can reach thermal equilibrium with the matter fields. 
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5.1 Single-hole Swiss-cheese 

As we have seen in Section |3.2|, both 2 <-> 2 and 1 2 processes can establish and maintain 



thermal equilibrium. Let us now apply the general conditions of Sections 3.2.1 and 3.2.2 
to our case. 

As we have already pointed out, scattering and annihilation processes involving strong 
interactions will establish thermal equilibrium between MSSM particles for temperatures 
T < a%g* ' Mp ~ 10 15 GeV. Let us now concentrate on the moduli. 

Small modulus $ 



From Section 4.1, we know that the largest coupling of the small canonical modulus <I> is 



with the non-abelian gauge bosons denoted by X: 



Vv i 

C<s>xx = g^xX^F^F^ ', g^xx ~ jj- ~ (5.1) 



Therefore according to (3.20), scattering or annihilation and pair production processes with 
two gravitational vertices like X + X <^ & + X + & <^ X + otX + X^X + X, can 
establish thermal equilibrium between <1> and X for temperatures: 

where Tj denotes the freeze-out temperature of the modulus. Taking the number of 



degrees of freedom g* to be 0(100), as in the MSSM, we find that fl|J) implies T > 5 x 10 8 
GeV for V ~ 10 15 , whereas T > 10 16 GeV for V ~ 10 4 . 12 In fact, for a typically large 
volume (V > 10 10 ) a more efficient 2 <-> 2 process is X + X <-» X + $ with one gravitational 



and one renormalisable vertex with coupling constant g. Indeed, according to ( 3.21 ), such 
scattering processes maintain thermal equilibrium for temperatures: 

1/2 

T > Tf ] = 9 * J^ P ~ 10 3 ^ for g* ~ 100 and 5 ~ 0.1 , (5.3) 
1 g z V V 

which for V ~ 10 15 gives T > 10 6 GeV while for V ~ 10 4 it gives T > 10 17 GeV. 

Finally, let us investigate the role played by decay and inverse decay processes of the 
form $ <-> X + V. We recall that such processes can, in principle, maintain thermal 
equilibrium only for temperatures: 

In V 

T > m$ ~ — M P , (5.4) 

because the energy of the gauge bosons is given by Ex ~ T and hence for T < m$ it 
is insufficient for the inverse decay process to occur. However, for T > m$ the process 
X+ X — > $ does take place and so one only needs to know the rate of the decay # — > V+ X 



in order to find out whether thermal equilibrium is achieved. According to (3.32) with 



"Recall that Mp here is the reduced Planck mass, which equals (87rGiv)" 1/2 = 2.4 x 10 18 GeV. 
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D ~ g^ xx /4Tr ~ V/4-7T, where we have also used ([5.1|), the condition for equilibrium is 
that: 

/ \ 1/3 / \ 1/3 

/ Vm$ \ / InV \ , . 

T <T eq = — m$ ~ -pr m$ = Km$ . (5.5) 

\Airgl /2 M P ) \W 7 / 

Hence thermal equilibrium between $ and X can be maintained by 1 <-> 2 processes only 
if re > 1 13 . However, estimating the total number of degrees of freedom as g* ~ 0(100), 
and writing the volume as V ~ 10 x , we obtain that k > 1 44> x > 55. Such a large value 
is unacceptable, as it makes the string scale too small to be compatible with observations. 
Therefore, we conclude that in LVS the small modulus $ never thermalises via decay and 
inverse decay processes. 

The final picture is the following: 

• For V of order 10 15 (10 10 ), as in typical LVS, from (|5.3| ) we deduce that the modulus 

is in thermal equilibrium with MSSM particles for temperatures T > ~ 10 6 

GeV (T > rj 2) ~ 10 11 GeV) due to X + X $ + X processes. 

• On the other hand, for V < 10 10 , as for LVS that allow gauge coupling unification, 
the main processes that maintain thermal equilibrium of the modulus with MSSM 
particles are purely gravitational: X +X <f>+X $+X or X +X ^ X +X 
and the freeze-out temperature is given by For example for V ~ 10 4 

M s ~ 10 16 GeV), is in thermal equilibrium for temperatures T > ~ 5 x 10 15 
GeV. 

We stress that this is the first example in the literature of a modulus that reaches thermal 
equilibrium with ordinary particles for temperatures significantly less than Mp, and so 
completely within the validity of the low energy effective theory. Note that we did not 
focus on the interactions of $ with other ordinary and supersymmetric particles, since the 
corresponding couplings, derived in Appendix [A|, are not large enough to establish thermal 
equilibrium. 

Finally, let us also note that, once the modulus $ drops out of thermal equilibrium, 
it will decay before its energy density can begin to dominate the energy density of the 
Universe, unlike traditional expectations in the literature. We will show this in more detail 



in Subsection S.3 



Large modulus x 



As summarised in Section 4.1, the coupling of the large modulus x with gauge bosons is 
given by 

£ x xx = 9 x xxxF^F^ v , g x xx ~ M p\ nV - ( 5 - 6 ) 



The exact value of k can be worked out via a more detailed calculation, very similar to the one that 



we will carry out in Section 6.4, It turns out that this value differs from the '~' estimate in (5.5) just by 
a multiplicative factor c 1//3 of 0(1). More precisely, c = 18(TT{r a })^ 3 ^' 2 e Kca ^ 2 Wo\ / 10gs and so, for natural 
values of all the parameters: Wo = 1, g a = 0.1, (r s ) = 5, K cs = 3, we obtain c 1 ^ 3 = 1.09. 
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Consequently purely gravitational 2^2 processes like X + X <-» x + X> X + x X + x, or 
X + X <-> X + X, could establish thermal equilibrium between x an d X for temperatures: 

T > Tj 1} = 5 y 6 M P (In V) 4/3 . (5.7) 

On the other hand, scattering processes like X + X <-» X + x with one gravitational and 
one renormalisable vertex with coupling constant g, could maintain thermal equilibrium 
for temperatures: 

1/2^ 

T > Tf ] = 9 * P (InV) 2 ~ 10 3 M P (In V) 2 , for ~ 100 and ff ~ 0.1 . (5.8) 
S 2 

Clearly, both t! and are greater than Mp and so we conclude that x can never 
thermalise via 2 <-> 2 processes. It is also immediate to notice that thermal equilibrium 
cannot be maintained by 1 <-> 2 processes, like x A + X, either. The reason is that, 



as derived in [17], for typical LARGE values of the volume V ~ 10 10 — 10 15 , the lifetime 
of the large modulus x is greater than the age of the Universe. Hence this modulus 
could contribute to dark matter and its decay to photons or electrons could be one of the 
smoking-gun signal of LVS. 



Furthermore, as can be seen from Section 4.1, the couplings of x to other MSSM 
particles are even weaker than its coupling to gauge bosons. So x cannot thermalise via 
any other kind of interaction. Finally, one can also verify that thermal equilibrium between 
X and $ can never be maintained via 1 <->• 2 and 2 «-> 2 processes involving only the moduli, 



which processes arise due to the moduli triple self-couplings computed in Appendix A. 3 
Therefore, x behaves as a typical modulus studied in the literature. 

5.2 Multiple-hole Swiss-cheese 



We shall now extend the results of Section 5.1 to the more general case of CY three- folds 



with one large cycle and several small ones. We shall not focus on explicit models since 
this is beyond the scope of our paper, but we will try to discuss qualitatively the generic 
behaviour of small moduli in the case of 'multiple-hole Swiss-cheese' CY manifolds. 



As we have seen in Section 4.1 , the couplings with MSSM particles of all the small cycles 



wrapped by MSSM branes have the same volume scaling as the corresponding couplings of 



the single small modulus in the CP^ 1 1 6 g , case. Moreover, in Section 5.1 we have learned 



that $ can thermalise via its interaction with gauge bosons. Hence, we conclude that the 



same arguments as in Section 5.1 can be applied for h\ \ > 2 and so all small cycles, that 



support MSSM chiral matter, reach thermal equilibrium with the gauge bosons. 



Note however that, as we already pointed out in Section 4.1, the situation may be 
more complicated in concrete phenomenological models due to the possibility that non- 
perturbative effects may be incompatible with MSSM branes, which are localized on the 
same 4-cycle [^3| . Whether or not such an incompatibility arises depends on the particular 
features of the model one considers, including the presence or absence of charged matter 
fields with non-vanishing VEVs. As a consequence of these subtleties, the issue of moduli 
thermalisation is highly dependent on the possible underlying brane set-ups. To gain 
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familiarity with the outcome, let us explore in more detail several brane set-ups in the 
case of only two small moduli. At the end we will comment on the generalization of these 
results to the case of arbitrary hix. 

We will focus on the case h± ; i = 3 with two small moduli t\ and T2, that give the 
volumes of the two rigid divisors Y\ and Ti- The results of Subsection [iTl] imply the 
following for the different brane set-ups below: 

1. If Ti is wrapped by an ED3 instanton and T2 is wrapped by MSSM branes: 

• n couples to MSSM gauge bosons with strength g ~ l/(VVMp) =4> t\ does not 
thermalise. 

• T2 couples to MSSM gauge bosons with strength g ~ \/V '/Mp =^ T2 thermalises. 

2. If Ti is wrapped by an ED3 instanton and Ti + T2 is wrapped by MSSM branes with 
chiral intersections on I^ 15 : 

• n couples to MSSM gauge bosons with strength g ~ \/V '/Mp =^ t\ thermalises. 

• T2 couples to MSSM gauge bosons with strength g ~ \/V/Mp =^ T2 thermalises. 

3. If Ti is supporting a pure SU(N) theory, that undergoes gaugino condensation, and 
Y2 is wrapped by MSSM branes: 

• t\ couples to MSSM gauge bosons with strength g ~ l/(VVMp) and to hid- 
den sector gauge bosons with strength g ~ W '/Mp =^ r% thermalises via its 
interaction with hidden sector gauge bosons. 

• T2 couples to MSSM gauge bosons with strength g ~ \fV /Mp and to hidden 
sector gauge bosons with strength g ~ l/(s/VMp) =^ T2 thermalises via its 
interaction with MSSM gauge bosons. 

Hence in this case there are two separate thermal baths: one contains n and the 
hidden sector gauge bosons at temperature Ti, whereas the other one is formed by 
T2 and the MSSM particles at temperature T2. Generically, we would expect that 
T\ 7^ T2 since the two thermal baths are not in contact with each other. 

4. If Ti is supporting a pure SU(N) theory, that undergoes gaugino condensation, and 
Ti + 1?2 is wrapped by MSSM branes with chiral intersections on IV 

• T\ couples both to MSSM and hidden sector gauge bosons with strength g ~ 
yV/Mp n thermalises. 



4 The coupling g ~ l/(y/VM P ) can be worked out by substituting the expression (14. 141) in (|A.7|). As 



pointed out in point 1 at the end of Subsection h. 2, the weakness of this coupling is due to the mixing term 



in (4.14) being highly suppressed by inverse powers of V. 

15 We assume that a single D7 brane is wrapping F2 in order to get chirality from the intersection with 
the MSSM branes. The same assumption applies throughout the paper everywhere we use the expression 
'chiral intersections on some divisor'. 
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• T2 couples to MSSM gauge bosons with strength g ~ W/Mp and to hidden 
sector gauge bosons with strength g ~ l/(s/VMp) =^ T2 thermalises via its 
interaction with MSSM gauge bosons. 

Unlike the previous case, now there is only one thermal bath, which contains both t\ 
and T2 together with the MSSM particles and the hidden sector gauge bosons, since 
in the present case t\ interacts strongly enough with the MSSM gauge bosons. 

We can now extend these results to the general case with h^i > 3 by noticing that a small 
4-cycle wrapped by MSSM branes will always thermalise via its interaction with MSSM 
gauge bosons. On the other hand, for a 4-cycle that is not wrapped by MSSM branes there 
are the following two options. If it is wrapped by an ED3 instanton, it will not thermalise. 
If instead it is supporting gaugino condensation, it will reach thermal equilibrium with the 
hidden sector gauge bosons. 

5.3 K3 Fibration 

Let us now turn to the issue of moduli thermalisation for K3 fibrations. As we have seen 



in Subsection [4.3| , there is an essential difference between the cases when the K3 fiber is 
stabilized at a large and at a small value. Let us consider separately each of these two 
situations. 

Large K3 fiber 



As we have already stressed in Subsection 4.3, in the case 'LV where the K3 divisor 
is stabilised large, the small modulus plays exactly the same role as the small modulus 
of the single-hole Swiss-cheese case, whereas both xi an d X2 behave as the single large 



modulus. Hence we can repeat the same analysis as in Subsection 5.1 and conclude that 
only will reach thermal equilibrium with the MSSM particles via its interaction with the 
gauge bosons. 

Small K3 fiber 

The study of moduli thermalisation in the case of small K3 fiber is more complicated. 
We shall first focus on CY three-folds with just one blow-up mode and later on will infer 
the general features of the situation with several blow-ups. 

K3 fibrations with h\ \ = 3 are characterised by two small moduli: t\ that gives the 
volume of the K3 divisor T\, and t s which is the volume of the rigid divisor T s . The 



canonically normalised fields \i an d $ are defined by ( 4.30| ) and ( [4.32| ). We recall that one 



has to be careful about the possible incompatibility of MSSM branes on T s with the non- 
perturbative effects that this cycle supports. Hence, to avoid dealing with such subtleties, 
below we will assume that the MSSM branes are not wrapping T s . Again, using the results 



of Subsection 5.1, we infer the following for the different brane set-ups below: 



1. If T s is wrapped by an ED3 instanton and T\ is wrapped by MSSM branes: 

• xi couples to MSSM gauge bosons with strength g ~ 1/Mp =^ xi does not 
thermalise. 
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• $ couples to MSSM gauge bosons more weakly than xi $ does not thermalise. 

2. If T s is wrapped by an ED3 instanton and T s + T± is wrapped by MSSM branes with 
chiral intersections on Ti'. 

• xi couples to MSSM gauge bosons with strength g ~ \/V/Mp =^ xi thermalises. 

• $ couples to MSSM gauge bosons with strength g ~ \/V jMp thermalises. 

3. If r ' s is supporting a pure S'JJ (iV) theory, that undergoes gaugino condensation, and 
Y\ is wrapped by MSSM branes: 



Xi couples to MSSM gauge bosons with strength g ~ 1/Mp and to hidden sector 
gauge boson with strength g ~ \fV jMp =^ xi thermalises via its interaction with 
hidden sector gauge bosons. 

• $ couples to MSSM gauge bosons more weakly than xi an d to hidden sector 
gauge bosons with strength g ~ \^V/Mp =^ $ thermalises via its interaction 
with hidden sector gauge bosons. 

In this case, two separate thermal baths are established: one contains xi> 3> an d the 
hidden sector gauge bosons at temperature Ti, whereas the other one is formed by 
the MSSM particles at temperature T2. Generically, we expect that T\ 7^ T2 since 
the two thermal baths are not in contact with each other. 

4. If r s is supporting a pure SU (N) theory, that undergoes gaugino condensation, and 
T s + Ti is wrapped by MSSM branes with chiral intersections on T\: 

• xi couples both to MSSM and hidden sector gauge bosons with strength g ~ 
W/Mp xi thermalises. 

• $ couples both to MSSM and hidden sector gauge bosons with strength g ~ 
VV/Mp =^ $ thermalises. 

Now only one thermal bath is established containing xi, 3>, the hidden sector gauge 
bosons and the MSSM particles, since both moduli interact with equal strength with 
the gauge bosons of the MSSM and of the hidden sector. 

It is interesting to notice that both moduli xi an d $ thermalise in all situations, except 
when the blow-up mode is wrapped by an ED3 instanton only. In this particular case, no 
modulus thermalises. It is trivial to generalise these conclusions for more than one blow-up 
mode and the MSSM still localised on the K3 fiber. 

On the other hand, if the MSSM is localised on one of the rigid divisors, then for the 
case of more than one blow-up mode one can repeat the same general conclusions as at the 



end of Subsection ^2, with in addition the fact that xi will always thermalise as soon as 
one of the blow-up modes thermalises. This is due to the leading order mixing between $ 
and any other small modulus, as can be seen explicitly in (4.33) and ( f4.34j ). 
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5.4 Modulini thermalisation 



The study of modulini thermalisation is straightforward since, as we have seen in Subsec- 



tion 4.4 , the canonical normalisation for the modulini takes exactly the same form as the 
canonical normalisation for the moduli. This implies that, after supersymmetrisation, the 
small modulino-gaugino-gauge boson coupling has the same strength as the small modulus- 
gauge boson-gauge boson coupling. Given that this is the relevant interaction for moduli 
thermalisation, we can repeat the same considerations as those in Subsections 5.1-^5 and 



conclude that the modulini thermalise every time, when their supersymmetric partners 
reach thermal equilibrium with the MSSM thermal bath. Note however that, if for the 
moduli the relevant processes are 2 <-> 2 interactions with gauge bosons, the crucial 2 *-* 2 
processes for the modulini are: 

• 2 «-> 2 processes with two gravitational vertices dominant for V < 10 10 : X + X «-> 

$ + 4>, x + x $ + $, x + $ x + $, x + $ <r+ x + x + x ^ x + x, 

X + X <-> X + X. 

• 2 *-* 2 processes with one gravitational and one renormalisable vertex dominant for 
V > 10 10 : X + ^ X + X, X + ^ X + X. 



6. Finite temperature corrections in LVS 

In this Section we study the finite temperature effective potential in LVS. We show that it 
has runaway behaviour at high T and compute the decompactification temperature T max . 
We also investigate the cosmological implications of the small modulus decay. By imposing 
that the temperature just after its decay (regardless of whether or not that decay leads to 
reheating) be less than T max , in order to avoid decompactification of the internal space, we 
find important restrictions on the range of values of the CY volume. 

6.1 Effective potential 

We shall now derive the explicit form of the finite temperature effective potential for LVS, 
following the analysis of moduli thermalisation performed in Section [5|. We will study in 
detail the behaviour of thermal corrections to the T = potential of the simple CPA 1169 j 
model, and then realise that the single-hole Swiss-cheese case already incorporates all the 
key properties of the general LVS. 

Single-hole Swiss-cheese 



As we have seen in Section 5.1, not only ordinary MSSM particles thermalise via Yang- 
Mills interactions but also the small modulus and modulino reach thermal equilibrium with 
matter via their interactions with the gauge bosons. Therefore, the general expression (^ 



for the 1-loop finite temperature effective potential, takes the following form: 

Vt- 1 °° P = (g B + ^f) + ^ (m% + m\ + Y. M s 2 OFT,i) + • (6-1) 
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We recall that (|6.l| ) is a high temperature expansion of the general 1-loop integral (|3.4j ), and 
so it is valid only for T 3> m$, m$, -Msoftv- The general moduli-dependent expression for 
the modulino mass-squared m| is given by ( 4,40 ) without the vacuum expectation value. 
On the other hand, in the limit 3> r s , m| can be estimated as follows: 

ml ~ TrM? - *£L^L „ (6 2 ) 



(6.3) 



For a s r s S> 1, the previous expression ( |6.2| ), at leading order, becomes: 



7T V \/2V 



It can be shown that the gaugino and scalar masses arising from gravity mediated SUSY 
breaking 16 are always parametrically smaller than m<j> and m^, and so we shall neglect 
them. Moreover we shall drop also the 0{T A ) term in ( |6.lD since it has no moduli depen- 
dence. Therefore, the relevant 1-loop finite-temperature effective potential reads: 

T 2 / 

rl-loop x I 2 , 2 



V T P = ^ (mi + m%) +... , (6.4) 
which using ( 4.40| ) and (6.3), takes the form: 



v i-ioo P _ T 2 f g s e K -Mf 



Wr 



+ ... , (6.5) 



Ai = 108A 2 s aj, X 2 = 3a 2 s A s W /V2. (6.6) 



24 V 7T 

with 



Given that the leading contribution in fl6.1| ), namely the C(T 4 ) term, does not bring in 
any moduli dependence, we need to go beyond the ideal gas approximation and consider 
the effect of 2-loop thermal corrections, as the latter could in principle compete with the 
terms in (^3J^). The high temperature expansion of the 2-loop contribution looks like: 

v 2-ioo P s _ t a (^ Kl g 2 MSSM + K 2 g% xx m\ + ^g% xx m\ + + ... , (6.7) 

where the re's are 0(1) coefficients and: 

• the 0{g\ ISSM ) contribution comes from two loops involving MSSM particles, 

• the 0{g% xx ) contribution is due to two loop diagrams with $ and two gauge bosons, 

• the 0(q\ xx ) contribution comes from two loops involving the modulino <£, the gaug- 
ino X and the gauge boson X, 



3 The contribution from anomaly mediation is subleading with respect to gravity mediation as shown in 
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all the other two loop diagrams give rise to subdominant contributions, and so 
they have been neglected. Such diagrams are the ones with $ or $ plus other 
MSSM particles, the self-interactions of the moduli and of the modulini, and two 
loops involving both $ and <3?. For example, the subleading contribution originating 
from the two-loop vacuum diagram due to the <& 3 self-interaction takes the form: 

dv 2-loops = T ^_h „ T 4 ^^. 
-f * mj V(lnV)' ! 



Note that in (6.7) we have neglected the 0(T 2 ) term since it is subleading compared to 
both the 0(T 4 ) 2-loop term and the 0(T 2 ) 1-loop one. Now, the relevant gauge couplings 
in (6/7), have the following moduli dependence: 



9mssm = 4tt/t s since we assume that the MSSM is built via magnetised D7 branes 
wrapping the small cycle. In the case of a supersymmetric SU(N C ) gauge theory with 
Nj matter multiplets, the coefficient K\ reads ]4(| : 



^ = ^{^-l)(N c + 3N f )>0. (6.8) 



9%xx ~ 9%xx ~ 1& as derived in (EH and O- 



Adding fl6.4| ) and (S.7) to the T = potential Vq, we obtain the full finite temperature 
effective potential: 

Vtot = V +T a (Kig 2 MSSM + K 2 glxx m l + K 34xx m l) + ^4 ( m * + m l) + - ■ ( 6 - 9 ) 

Despite the thermalisation of $ and 3>, which in principle leads to a modification of Vtot 
compared to previous expectations in the literature, we shall now show that the thermal 
corrections due to $ and $ are, in fact, negligible compared to the other contributions in 
fl6.9p , everywhere in the moduli space of these models. In particular, the 2-loop MSSM 
effects dominate the temperature-dependent term. 17 

Let us start by arguing that the 0(T 4 ) corrections arising from the modulus and 
the modulino <I> are subleading compared to the 1-loop 0(T 2 ) term. Indeed, the relevant 
part of the effective scalar potential (|6.9|) may be rewritten as: 



T (^g^xx™® + K 3g$xx m H>) ~ T ( m $ + m i.) T T7T ' ( 6 - 10 ) 




where the « inequality is due to the fact that our effective field theory treatment makes 
sense only at energies lower than the string scale M s . Therefore, we can neglect the effect 
of 2-loop thermal corrections involving $ and So we see that, although the interactions 
of $ and $ with gauge bosons and gauginos are strong enough to make them thermalise, 



17 Note that this is consistent with the results of in the context of the O'KKLT model, where it was 
also found that the T-dependent contribution of moduli, that were assumed to be in thermal equilibrium, 
is negligible compared to the dominant contribution of the rest of the effective potential. 
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they are not sufficient to produce thermal corrections large enough to affect the form of the 
total effective potential. Let us also stress that this result is valid everywhere in moduli 
space, i.e. for each value of m\ and m|, not just in the region around the zero-temperature 
minimum. 

We now turn to the study of the general behaviour of the 1-loop OiT 2 ) term aris- 
ing from and <I>. We shall show that it is always subdominant compared to the zero- 
temperature potential ( |2.17 ), and so it can be safely neglected. In fact, the two relevant 
terms (E17D and Q can be written as (ignoring the subleading loop corrections in Vq): 



T 2 ( 

with 



2 i 2 
771$ + 771; 



K c 



M p 



8tt 



pi = 36a 4 s A 2 s , 



P2 



-2a„T. 



V 



4a s A s W 



p 2 A : 



V 2 



P3 = W 2 /6, 



(6.11) 



(6.12) 



and 



2y/2 T 2 V^T S 



3a 2 



+ 



M 2 



A 2 = l + 



T 2 V, 



4^/2 



Mp 



A 3 = ^ + 



T 2 V 



( «L ) 



<1 



<1 




where the appearance of the Kaluza-Klein scale comes from the assumption that the MSSM 
branes are wrapping the small cycle t s : 

M s M P 



M K k 



1/4 



1/4 



(6.13) 



Therefore, we can see that the 1-loop OiT 2 ) thermal corrections can never compete with 
Vq for temperatures below the compactification scale Mkk < M s , where our low energy 
effective field theory is trustworthy. Once again, we stress that the previous considerations 
are valid in all the moduli space (within our large volume approximations) and not just 
in the vicinity of the T = minimum. We have seen that the only finite-temperature 
contribution that can compete with Vq is the 2- loop T 4 gj^ SSM term, and so we can only 
consider from now on the following potential: 



Vtot = Vq + 4ttki- 



Kc 



8tt 



A, 



'r,e 



-2a»T s 



V 



V 2 



v 



4^£i (TV 
~r~T \Mp~J 



M P , 



(6.14) 



valid for temperatures T 3> Msoft, and with the constants given in ( p,18| ) and (|6.S]) 18 . 
We realize that the leading moduli-dependent finite temperature contribution to the effec- 
tive potential comes from 2-loops instead of 1-loop. This, however, does not mean that 
perturbation theory breaks down, since 1-loop effects still dominate when one takes into 
account the moduli independent 0(T 4 ) piece that we dropped. 

Now, from ( 6.14| ) it is clear that the thermal correction cannot induce any new In- 
dependent extremum of the effective potential. Its presence only leads to destabilization 



3 For convenience, here we have redefined ki = 8ivKig 3 
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of the T = minimum at a certain temperature, above which the potential has a runaway 
behaviour. Therefore, we are led to the following qualitative picture. Let us assume that at 
the end of inflation the system is sitting at the T = minimum. Then, after reheating the 
MSSM particles thermalise and the thermal correction T 4 gj^ SSM ~ T 4 /t s gets switched 
on. As a result, the system starts running away along the r s direction only, since Vt does 
not depend on V. However, as soon as r s becomes significantly larger than its T = VEV, 
the two exponential terms in (6.14) become very suppressed with respect to the 0(V~ 3 ) 
a' correction (the v term). Hence, the potential develops a run-away behaviour also along 
the V-direction, thus allowing the Kahler moduli to remain within the Kahler cone. 



In Section p.2| , we shall compute the decompactification temperature, at which the 
T = minimum gets destabilised. Hence we shall focus on the region in the vicinity of the 
zero-temperature minimum, where the regime of validity of the expression flOg ) takes the 
form: 

MsOFT < T < MkK T77-TT < — < ^ ma - (6-15) 



VlnV M P VVr s 1/4 ' 
10 14 allows low energy SUSY, we get Msoft 



In the typical LVS where V ~ 10 i4 allows low energy SUSY, we get Msoft ~ 10 3 GeV 
and Mkk ~ 10 11 GeV; thus, in that case, eq. (6.14) makes sense only for energies 10 3 



GeV « T « 10 11 GeV. On the other hand, for LVS that allow GUT string scenarios, 
V ~ 10 4 , which implies M SO ft ~ 10 13 GeV and M KK ~ 10 16 GeV; thus, in that case, 
flOp becomes 10 13 GeV <T« 10 16 GeV. 

General LARGE Volume Scenario 



As we have seen in Section 2.3, one of the conditions on an arbitrary Calabi-Yau to 
obtain LVS, is the presence of a blow-up mode resolving a point-like singularity (del Pezzo 
4-cycle). The moduli scaling of the scalar potential, at leading order and in the presence 
of N srna ii blow-up modes r s< , i = 1, N sma a, is still of the form ( 2.17] ) (neglecting loop 
corrections) : 



Vh 



9 s e 



K c 



M P 



8tt 



-^V small 

£ 

. i=i 



A. 



"To. e 



-a s .T s , 



V 



V 2 



+ 



V 3 



(6.16) 



All the other moduli which are neither the overall volume nor a blow-up mode will appear 
in the scalar potential at subleading order. Moreover, due to the topological nature of r Sj j, 
VJfT: Vi = l,...,N small pi. 



As derived in Section 4.1, these blow-up modes correspond to the heaviest moduli and 



modulini, which play the same role as $ and $ in the single-hole Swiss-cheese case. Hence 
the leading order behaviour of the mass-squareds of the blow-up moduli r s . and the cor- 
responding modulini f Sj are still given by (|6.3| ) and ( [4.40 ) \fi = 1, N sma u. Therefore we 
can repeat the same considerations made in the previous paragraph and conclude that, for 
a general LVS, the 1-loop OiT 2 ) thermal corrections are always subdominant with respect 
to Vo for temperatures below the compactification scale 19 . The only finite-temperature 

19 As we have seen in Section |5.2[ if all the T Si are wrapped by ED3 instantons then they do not thermalise. 
Only the moduli corresponding to 4-cycles wrapped by MSSM branes would then thermalise but, since they 
are lighter than the ED3 moduli, our argument is still valid. The same is true for all the possible scenarios 



outlined for the K3 fibration case in Section 5.3 
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contribution that can compete with Vq is again the 2-loop 1 Smssm term. 



6.2 Decompactification temperature 

As we saw in the previous subsection, the finite temperature corrections destabilize the large 
volume minimum of a general LVS. In this subsection we will derive the decompactification 
temperature T max , that is the temperature above which the full effective potential has no 
other minima than the one at infinity. 

Before performing cL mor6 precisB calculation of T max: let us present a qualitative 
argument that gives a good intuition for its magnitude. Let us denote by Vb the height of 
the potential barrier that separates the supersymmetric minimum at infinity from the zero 
temperature SUSY breaking one. Now, in order for the moduli to overcome the potential 
barrier and run away to infinity, one needs to supply energy of at least the same order of 
magnitude as Vb- In our case, the source of energy is provided by the finite-temperature 
effects, which give a contribution to the scalar potential of the order Vp ~ T A . Hence a 
very good estimate for the decompactification temperature is given by T max ~ V^ 4 . 

It is instructive to compare the implications of this estimate for the KKLT and LVS 
cases. In the simplest KKLT models the potential reads: 

Vkklt = Xi AalVo^, (6.17) 

where Ai and A2 are constants of order unity. The minimum is achieved by fine tuning the 
flux parameter Wq ~ re _ar and so the height of the barrier is given by 

W 2 

V b ~ (V KKL t} ~ ~^M P ~ ml /2 Mj,, (6.18) 

where we have used the fact that V = t 3 / 2 and m 3/ / 2 = WqMp/V. Therefore the de- 
compactification temperature becomes T max ~ ^Jm^Mp ~ 10 10 GeV, as estimated in 

@- 

In the case of LVS, the height of the barrier is lower and so we expect a lower decom- 
pactification temperature T max . Indeed, to leading order the potential is given by 

e -2a s T 3 e -a 3 T s W2 

V L vs = XiV^s— X 2 W t s ^- + A 3 ^f- (6.19) 

with Ai, A2 and A3 being constants of order one, as reviewed in Section ||[ The minimum 
is achieved for natural values of the flux parameter Wq ~ 0(1) and at exponentially large 
values of the overall volume V ~ WQ^/rle aaTs . Hence the height of the barrier can be 
estimated as: 

V b ~ (Vlvs) ~ y§-M£ ~ rn\ j2 M P , (6.20) 
which gives a decompactification temperature of the order: 



m | /2 Mp)""~^. (6.21) 



1/4 M P 

Let us now turn to a more precise computation. Without loss of generality, we shall 
focus here on the effective potential ( |6.14| ), valid for the single-hole Swiss-cheese case, and 
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look for its extrema. Given that the thermal contribution does not depend on the volume, 
the derivative of the potential with respect to V gives the same result as in the T = case: 



OVtot 
dV 



V* = ^A(r s )^e a 



where 20 



A(r s ) = 1 



(Ts) 



3/2 



(6.22) 



(6.23) 



and (t s ) ~ (AXu/fi 2 ) 2/3 is the T = VEV of t s . Substituting ([6.22j ) in the derivative of 
Vtot with respect to r s and working in the limit a s r s 3> 1, in which one can neglect higher 
order instanton corrections, we obtain: 



dVroT 







4"7TK 



fie 



3a 3 T a 



V=V, 



A 2 a., 



Mp 



A(t s ) 2 + 2A(t s ) -1 = 0. 



(6.24) 



Notice that at zero temperature (6.24) simplifies to A(t s ) = 1/2, which from ( 6.23| ) correctly 
implies r s = (t s ). Now, since equation ( |6.24 ) is transcendental, one cannot write down an 
analytical solution, that gives the general relation between the location of the t s extrema 
and the temperature. Nevertheless, we will see shortly that it is actually possible to extract 
an analytic estimate for the decompactification temperature. To understand why, let us 
gain insight into the behaviour of the function on the LHS of ( |6.24 ) by plotting it and 
looking at its intersections with the r s -axis. 



We plot the LHS of equation ( |6.24 ) on Figure [5] for several values of the temperature; 



T increases from right to left. From this figure it is easy to see that the temperature- 
dependent correction to Vtot behaves effectively as an up-lifting term. Namely, the finite- 
temperature contribution lifts the potential, giving rise to a local maximum (the right 
intersection with the t s axis) in addition to the T = minimum (the left intersection). As 
the temperature increases, the maximum increases as well and shifts towards smaller values 
of t s . On the other hand, the minimum remains very close to the zero-temperature one at all 
temperatures. Clearly, the decompactification temperature T max is reached when the two 
extrema coincide. The key observation here is that this happens in a small neighborhood 
of the T = minimum, located at (t s ) ~ (4Az;/^ 2 ) 2 ^ 3 . 

In view of the considerations of the previous paragraph, to find an analytic estimate 
for T max we shall utilize the following strategy. We will Taylor-expand the function F(r s ), 
defined by the LHS of equation ( |6.24 ), to second order in a small neighborhood of the point 



t s = (t s ). Then we will use the resulting quadratic function f(S), where 5 = t s — (r s ), 
as an approximation of F{r s ) in a larger neighborhood and will look for the zeros of 
f(6). Requiring that the two roots of f(5) coincide, will give us an estimate for the 
decompactification temperature. Clearly, this procedure is not exact. In particular, the 



20 We discard the solution with the positive sign in front of the square root in ( |6.23| ) since, upon its 
substitution one finds that the other extremum condition, 8Vtot / dr s = 0, does not have any solution. 
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dV T0 - 




Figure 5: The LHS of eq. ( |6.24 ) is plotted versus r s . The temperature increases from right to 
left. The straight line represents the zero temperature case. The other values of the temperature 
are T/M P = 0.8 • KT 10 , 1.0 • lO" 10 , 1.2 • 10~ 10 , 1.4 • lO" 10 . To obtain the plots we used the 
following numerical values: £ = 1.31, A s = 1, Wo = 1, a s = tt/4, e Kcs = 8iT/g s , g s — 0.1, N c = 5, 
Nf = 7. With these values one has that (r s ) = 41.55 and (V) = 7.02 • 10 13 , which implies that 
T max = 1.58 • 10" 10 A/ P ~ 3.79 • 10 8 GeV according to (|6.30[) . Note that the numerically found value 
of the decompactification temperature is T max ^ num = 1.20 • 10 _10 A/p. 



function F(t s ) is better approximated by keeping higher orders in the Taylor expansion. 
In our case, we have checked numerically that a really good approximation is obtained by 
going to at least sixth order. However, in doing so one again ends up with an equation that 
cannot be solved analytically. So the key point is that the systematic error introduced by 
the quadratic approximation is rather small (we have checked that the analytical results 
obtained by following the above procedure are in very good agreement with the exact 
numerical values). 

Now let us substitute t s = (t s ) + 5 in ( 6.24 ) and read off the terms up to order 5 2 . 
The result is: 

a5 2 + b5 + c = 0, (6.25) 



where the corresponding coefficients, in the limit a s {r s ) 3> 1, take the form: 

f a ~ \7a 2 s + i|iA 2 a s , 

I b ~ 3Ta s - 9X 2 a s (T s ), c~T, 



(6.26) 



and we have set ^ 

T = 4kk 1 (-^-\ fxe 3a ^K (6.27) 

Finally, to find the decompactification temperature, we require that the two solutions 5± 
and 62 coincide: 

5 1 = 5 2 b 2 -4ac = 0, (6.28) 
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which, for a s (r s ) 3> 1, gives: 

T max = 3(V2-l)X"(r s ) <=> ^ = ?!^ffi e -*.W4 (6.29) 
Notice that we can rewrite the decompactification temperature in terms of V as: 

T ^ = 32 ir A^H^~ M P => T max ^ [m 3/2 M P ) ~ ^ , (6.30) 

where we have used the relation between the T = VEV of the volume and (t s ), which 
is given by ( 6.22; ) with t s = (t s ) and A = 1/2. It is reassuring that ( 6.30 ) is of the same 



form as the result ( 6.21 ), obtained from the intuitive arguments based on the height of the 
potential barrier. 

6.3 Small moduli cosmology 



Clearly, the decompactification temperature ( |6.30 ) sets an upper bound on the temperature 
in the early Universe, in particular on the reheating temperature, T^ H , at the end of 
inflation. We will investigate now how this constraint affects the moduli thermalisation 

21 



picture studied in Subsection 5.1 



Recall that there we derived the following: 

• For small values of the volume (V < 10 10 ), the freeze-out temperature for the small 
modulus $ is given by ([5^ ): T§ v ~ M P V~ 2/3 . 

• For large values of the volume (V > 10 10 ), the freeze-out temperature for <I> is given 
by (pD: Tj v ~ 10 3 MpV- 1 . 

Note also that, in both cases, the condition Tf < < T max has to be satisfied in order 
for the modulus to reach equilibrium with the MSSM thermal bath. Now, for small values 
of V we have that: 

T V 2 / 3 

- L max v 



T SV V 3/4 



V~ 1/u < 1, (6.31) 



which implies that $ actually never thermalises. On the other hand, for large values of V 
we have that (writing V ~ W x ): 

T m ax V / 



Tf v 10 3 



10 a>/4-3 >1 x> n. (6.32) 



Hence, for V > 10 12 , $ can reach thermal equilibrium with the MSSM plasma, as long as 
Tft H is such that Tj V < T^ H < T max . Let us stress, however, that if Tp[ H < Tj V the 
modulus will never thermalise even though Tf v < Tmax- Note that, since the temperature 
Tfi H depends on the concrete realization of inflation and the details of the initial reheating 
process, its determination is beyond the scope of the present paper. So we will treat it as 
a free parameter, satisfying only the constraint T^ H < T max . 

We would like now to study the cosmological history of $ which, in our case, presents 
two possibilities: 

21 Similar considerations apply for the more general multiple- hole Swiss-cheese and K3 fibration cases. 
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1. The modulus $ decays at the end of inflation being the main responsible for initial 
reheating. We may envisage two physically different situations where this could 
happen: in one case, $ is the inflaton and it decays at the end of inflation. In 
the other case, <1> is not the inflaton, but it starts oscillating around its VEV when 
the inflaton is still driving inflation by rolling down its flat potential. In this case, 
the decay of <3? occurs just after the slow-roll conditions stop being satisfied and the 
inflaton reaches its VEV. 

After $ decays, its energy density is converted into radiation. The decay products 
thermalise rapidly and re- heat the Universe to a temperature Trh = T^ H . The 
latter can be computed by noticing that the <3? energy density p$ ~ T\_ iXX M^ > will 
be converted into radiation energy density pp ~ g*T 4 . Hence Tp\ H can be obtained 
by comparing T^^xx with the value of H, given by the Priedmann equation for 
radiation dominance: 

^ XX 16vr M P - U 9 * M P 

n / InV \ 1/2 (In V) 3/2 M P 

Trh — ( — — J m$ = ^—>-f. 6.33 

In order for this picture to be compatible with the presence of a decompactification 



temperature ( 6.30 ), that sets the maximal temperature of the Universe, we need to 



require that T^ H < T max . As we shall see in Subsection p-4 this requirement can be 



translated into a constraint on the values that the internal volume can take. 

2. The modulus $ is not the main source of initial reheating, which we suppose to be the 
inflaton. After the inflaton decays, the Universe is re-heated to a temperature T^ H 
and an epoch of radiation dominance begins. The modulus $ will only thermalise if 
V > 10 12 and Tf v < T^ H . However, Tf v is rather close to 

T ma x and so, even when 

$ thermalises, it will drop out of equilibrium very quickly at T^ v . Then, for general 
values of V, the modulus <3? will decay out of equilibrium at a temperature Tp> < T^ H - 
As we shall show below, this decay will occur during radiation domination, since 
Td > Tdom, with Tdom being the temperature at which the modulus energy density 
would dominate over the radiation energy density. So the temperature Tp> at which 
$ decays, is still given by ( |6.33| ) upon replacing Tp[ H with Tp>: 

(In vf 2 M P 

Td = <6 ' 34) 

Note that the above expression satisfies Tp> < T^ V ' LV , as should be the case for 
consistency. Another important observation is that ( |6.34| ) is also the usual expression 
for the temperature Tr#, to which the Universe is re- heated by the decay of a particle 
releasing its energy to the thermal bath. In other words, for us Tpn = Tp> since 
the modulus $ decays during radiation domination. On the contrary, if a modulus 
decays when its energy density is dominating the energy density of the Universe, then 



-43- 



Td < Trh and the decay produces an increase in the entropy density S, which is 
determined by: 

A "ff~(^) 3 - (6 - 35) 
As already mentioned, since for us Trh = Tjj, the decay of $ does not actually lead 
to reheating or, equivalently, to an increase in the entropy density, given that from 
( ggD we have A = 1. As a consequence, <I> cannot dilute any unwanted relics, like for 
example the large modulus x which suffers from the cosmological moduli problem. 22 

To recapitulate: in the present case 2, we have the following system of inequalities: 

for V < 10 12 : T dom < T D < 1% H < T max , (6.36) 
for V > 10 12 : T dom < T D < Tf V < T° RH < T max . (6.37) 

As in case 1 above, the condition Tjy < T max implies a constraint on V, that we will 
derive in Subsection 6A. We underline again that this condition is necessary but 
not sufficient, since for us T RH is an undetermined parameter. In concrete models, 
in which one could compute T RH , the condition T RH < T max might lead to further 
restrictions. 

Let us now prove our claim above that, when the modulus $ is not responsible for the 
initial reheating (case 2), it will decay before its energy density begins to dominate the 
energy density of the Universe. will start oscillating around its VEV when H ~ m$ at 
a temperature T osc given by: 

T osc ~ g* 1/4 y/m$M P . (6.38) 

The energy density p$, stored by <3?, and the ratio between p§ and the radiation energy 
density at T osc read as follows: 



Pr) 



ml(T s ) 2 (r s ) 2 
~ TIT- ~ -Tp- ( 6 - 39 ) 

T OS c y* ± osc 1V1 P 



By definition, the temperature Tdom-, at which />$ becomes comparable to p r and hence $ 
begins to dominate the energy density of the Universe, is such that: 



Pt> 
Pr 



~ 1. (6.40) 



Now, given that p§ redshifts as T 3 whereas p r scales as T 4 , we can relate T^om with T c 



- dom 





T ( P*\ 




~ J- osc 1 1 




T dom \Pr J 



T dom ~ 9* 1/4 \^-V^*M P . (6.41) 

rp iVl r> 

■i- osc JT 



We shall show now that T^ om < To with Tjj being the decay temperature during radiation 
dominance, which is obtained by comparing H with T^^xx'- 

T D ~ ff ; 1/4 \AWMp. (6.42) 



22 This kind of solution of the cosmological moduli problem, i.e. dilution via saxion or modulus decay, is 
used both in [43) and in Q. 



-44- 



The ratio of (|6.42[) and |T4l]) gives 



7d om y 7 ^" (r s ) 2 ' 
Using that T^ X x ~ Vm|Mp 2 and (t s ) ~ 10M S ~ 10M P V~ 1/2 , the last relation becomes: 



T D _.(lnV)^V >1 for V>1Q , 5 (644) 



T dom 100 

Hence, we conclude that To > T^ om and, therefore, $ decays before it can begin to dom- 
inate the energy density of the Universe. The main consequence of this is that <5 cannot 
dilute unwanted relics via its decay. 

6.4 Lower bound on V 

As we saw in the previous Subsection, there are two possible scenarios for the cosmological 
evolution of the small modulus <J>. However, since the RHS of ( |6.33|) and ( |6.34| ) coincide, 
in both cases the crucial quantity is the same, although with a different physical meaning. 
Let us denote this quantity by T* ~ (r$Mp) 1//2 . We shall impose that T* < T max and 
shall show below that from this requirement one can derive a lower bound on the possible 
values of V in a general LVS. Before we begin, let us first recall that: 



1. If $ is responsible for the initial reheating via its decay, then T* = TR 



RH ' 



2. If decays after the original reheating in a radiation dominated era, then T* = Tjj < 

RH' 

Regardless of which of these two situations we consider, T* is the temperature of the 
Universe after $ decays. Then, in order to prevent decompactification of the internal 
space, we need to impose T* < T max . In general, this condition is necessary but not 
sufficient because in case 2 one must ensure also that T^ H < T max . This is a constraint 
that we cannot address given that in this case T^ H is an undetermined parameter for us. 
Let us now compute T* precisely. We start by using the exact form of the decay rate 

2 3 

T ^ XX ~ -M^Mf' (6 - 45) 

where 

25/4 /q 

9*xx = T^u^V- (6-46) 

The mass of $ is given by: 

m, = VP 2 -^f^M P , (6.47) 

where we are denoting with P the prefactor of the scalar potential: P = g s e Kcs /(8ir). From 
the minimisation of the scalar potential we have that 

a s (r s ) = In (pV) = lnp + In V, (6.48) 
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where 



and so 



P 



l2y/2a a A a 



Oil) 



a s (T s ) ~ In V, 



Therefore, the decay rate T^^xx turns out to be: 

•3 



rv 



p3/2 3W$(lnVy M P 
V2tt(t s }^ 2 V 2 ' 



(6.49) 
(6.50) 

(6.51) 



Finally, in order to obtain the total decay rate, we need to multiply T^xx by the total 
number of gauge bosons for the MSSM Nx = 12: 



v tot p3 /2 36W 3 (In Vf Mp 
^ XX ^vr(r s )3/2 V^- 



(6.52) 



Now, we can find T* by setting 4 {T% ( ^ X x) 2 equal to 3H 2 , with H read off from the 
Friedmann equation for radiation dominance: 



r ,.(4^) ,/4 ^^.p 3 /^(») 1/4 (W!*. ( , 53) 

\n 2 g*J v vr \g*J (r s ) 3 / 4 V 

We are finally ready to explore the constraint T* < T max . Recall that the maximal tem- 
perature is given by the decompactification temperature ( |6.30| ): 



(v^-l) 
4\/2 



1/4 



V 3 / 4 



M P . 



(6.54) 



Let us now consider the ratio T max /T* and impose that it is larger than unity (using 
g*(MSSM) = 228.75): 



R 



V 1 / 4 



T* (InV) 3 / 2 
where we have defined: 



with 



J 



(V2-l)g» 
80^2 



1/4 



11/ 



<r.> 



11/8 



6W 



2W 



-1/4 



8.42 

^V 4 



e -X cs /2 for 



9s 



0.1, 



(6.55) 



(6.56) 



and in fl6.55| ) we have set J = 1. In fact, from ( |6.8j ), we find that in the case of SQCD with 
N c = 3 and Nf = 6, ki = 2.625. However for the MSSM we expect a larger value of Ki 
which we assume to be of the order k% = 10. Then for natural values of K cs like K cs = 3 23 , 



23 The dependence of the Kahler potential on the complex structure moduli can be worked out by com- 
puting the different periods of the CY three-fold under consideration. As derived in [Q, for the simplest 
example of a CY manifold with just one complex structure modulus U (the mirror of the quintic), for 
natural values of U, \U\ ~ 0(1) K CB ~ 0(1). 
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15 



Figure 6: Plots of the ratio R = T max /T* as a function of V = 10 x and the parameter c rn m < c < 
Cmax as defined in (3.55), (6.57) and (6.58|). In the left plot, the red surface is the constant function 



R = 1, whereas in the right plot the black line denotes the curve in the (x,c)-plane for which R = 1. 




Figure 7: Plot of the R = 1 curve in the (a;,c)-plane. The shaded region represents the phe- 
nomenologically forbidden area, in which the values of x and c are such that R < 1 O- T max < T» . 



from ( 6.56j ), we find J = 1.05. Let us consider now the maximum and minimum values 
that the parameter c can take for natural values of (t s ) and Wq: 



-- 100 
0.01 



w 0)max = 100 



10 4 



10" 



(6.57) 
(6.58) 



Now writing V ~ 10 x , R becomes a function of x and c. Finally, we can make a 3D 
plot of R with c m i n < c < c max and 2 < x < 15, and see in which region R > 1. This 
is done in Figure ||. In order to understand better what values of V are disfavoured, we 
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also plot in Figure 0, as the shaded region, the region in the (x,c)-plane below the curve 
R = 1, which represents the phenomenologically forbidden area for which T max < We 
conclude that small values of the volume, which would allow the standard picture of gauge 
coupling unification and GUT theories, are disfavoured compared to larger values of V, that 
naturally lead to TeV-scale SUSY and are thus desirable to solve the hierarchy problem. 
In Table ||, we show explicitly how the lower bound on the volume, for some benchmark 
scenarios, favours LVS with larger values of V. 

From the definition ( |6.55 ) of the parameter c, 





R > 1 44> T max > T* 


c = 


4 


Vx 


c = 


3 


x > 2.1 


c = 


2 


x > 3.8 


c = 


1 


x > 5.9 


c = 


0.5 


x > 7.6 


c = 


0.1 


x > 11.3 


c = 


0.05 


x > 12.8 


c = 


0.01 


x > 16.1 



Table 4: Lower bounds on the vol- 
ume in the string frame V s ~ 10 x ~ 3 / 2 
for some benchmark scenarios. 



Hence if we write Ve ~ 10 : 
General LARGE Volume Scenario 



it is interesting to notice that for values of (t s ) far 
from the edge of consistency of the supergravity ap- 
proximation (t s ) ~ O(10), c should be fairly large, 
and hence the bound very weak, for natural values 
of Wo ~ 0(1), while c should get smaller for larger 
values of Wo that lead to a stronger bound. In addi- 
tion, it is reassuring to notice that for typical values 
of V ~ 10 15 , 

T m ax > except for a tiny portion of 
the (x,c)-space. It also important to recall that the 
physical value of the volume as seen by the string is 
the one expressed in the string frame V s , while we are 
working in the Einstein frame where V s 
then we have that V s = 10 x ^ 3 ^ 2 , upon setting g s = 0.1. 



gl /2 V E . 



Let us now generalise our lower bound on V to the four cases studied in Subsections 5.2 



and 5.3 for the multiple-hole Swiss-cheese and K3 fibration case (focusing on the small K3 
fiber scenario) respectively. 

First of all, we note that, since in all the cases the 4-cycle supporting the MSSM is 
stabilised by string loop corrections ||, we can estimate the actual height of the barrier 
seen by this modulus as (see (|2.22D ): 



Vb 



wl_ 



(6.59) 



where we are generically denoting any small cycle (either a blow-up or a K3 fiber divisor) 
as r, given that the values of the VEV of all these 4-cycles will have the same order of 
magnitude. Then setting ~ T max /r, we obtain: 



T ~ 

^ max 



v 3 ' 



We notice that ( 6.60 ) is a bit lower than ( 6.30 ) but the two expressions for T„ 
same leading order V-dependence. 



(6.60) 



share the 



Let us now examine the 4 cases of Subsection 5^ in more detail, keeping the same 
notation as in that Subsection, and denoting as $ the small modulus of the single-hole 
Swiss-cheese scenario studied above: 
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1. The relevant decay is the one of T2 to MSSM gauge bosons. The order of magnitude 
of the mass of T2 is: 

m T2 V 2 r 2 » ( 6 - 61 ) 

and so T2 is lighter than and, in turn, T* will be smaller. In fact, plugging (|6.6l| ) 
in ( 16.45 ) , we end up with (ignoring numerical prefactors): 



(lnV) 3 / 2 ^ 3/2 , , 

* - y r9 /4° • ( 6 - 62 ) 



Hence we obtain 



T V 1 / 4 -7- 19 / 8 
R{1) _ Irnax = (1 ) V (1) ^ T 

^ (lnV) 3/2 ^0 V 7 



Comparing this result with ( 6.55| ), we realise that i^ 1 ) ~ i?r and so the lower bound 



on V turns out to be less stringent. The final results can still be read from Table [I] 
upon replacing c with . 

2. The relevant decay is the one of n to MSSM gauge bosons since m Tl ~ m$, and so 
T\ is heavier than ti. Therefore will still be given by ( |6.53| ). Hence we obtain 

T V 1 / 4 t 7 / 8 

R<$ = ±™!L = cW - , l0 with C ( 2 )~^. (6.64) 
T* (lnV) 3/2 W V ; 

Comparing this result with ( |6.55| ), we realise that ~ Rt~ 1 / 2 and so the lower 
bound on V turns out to be more stringent. The final results can still be read from 
Table |] upon replacing c with . 

3. The relevant decay is the one of t\ to hidden sector gauge bosons. Hence we point 
out that the considerations of case 2 apply also for this case. 

4. The relevant decay is the one of t\ to MSSM gauge bosons, and so we can repeat the 
same considerations of case 2. 

The final picture is that for all cases the V-dependence of the ratio T max /T* is the same as 
in ( 6.55| ). The only difference is a rescaling of the parameter c. Thus we conclude that, as 



far as the lower bound on V is concerned, the single-hole Swiss-cheese case shows all the 
qualitative features of a general LVS. 

Finally, we mention that in the case of a K3 fibration with small K3 fiber, cases 2, 



3, and 4 of Subsection 5.3 have the same behaviour as case 2 of the multiple-hole Swiss- 
cheese, so giving a more stringent lower bound on V. We should note though that this 
lower bound does not apply to case 1 of Subsection |5,3|, since both of the moduli have 
an Mp-suppressed, instead of M s -suppressed, coupling to MSSM gauge bosons. However, 
these kinds of models tend to prefer larger values of V (due to the fact that a s = 2tt for an 
ED3 instanton) which are not affected by the lower bound that we derived. 
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7. Discussion 



Let us now discuss some of the possible applications of these results, as well as directions 
for future work. As we have emphasized throughout the paper, there are two kinds of LVS, 
depending on the magnitude of the value of the internal volume V. Their main cosmological 
characteristics are the following: 

LV case 

In this case the volume is stabilised at large values of the order V ~ 10 15 which allows 
to solve the hierarchy problem yielding TeV scale SUSY naturally. Here are the main 
cosmological features of these scenarios: 

• The moduli spectrum includes a light field \ related to the overall volume. This field 
is a source for the cosmological moduli problem (CMP) as long as M s < 10 13 GeV, 
corresponding to V > 10 10 . In fact, in this case the modulus \ is lighter than 10 
TeV and, coupling with gravitational strength interactions, it would overdose the 
Universe or decay so late to ruin Big-Bang nucleosynthesis. There are two main 
possible solutions to this CMP: 

1. The light modulus x gets diluted due to an increase in the entropy that occurs 
when a short-lived modulus decays out of equilibrium and while dominating the 
energy density of the Universe p3| . (b"f |; 

2. The volume modulus gets diluted due to a late period of low energy inflation 



caused by thermal effects [ 45 1 . 



Assuming this problem is solved, the volume modulus becomes a dark matter candi- 
date (with a mass m ~ 1 MeV, if V ~ 10 15 ) and its decay to e + e _ could be one of 
the sources that contribute to the observed 511 KeV line, coming from the centre of 
our galaxy. 24 The light modulus x can & lso decay into photons, producing a clean 
monochromatic line that would represent a clear astrophysical smoking-gun signal 



for these scenarios [17]. We point out that in the case of K3 fibrations, where the K3 
fiber is stabilised large ||, the spectrum of moduli fields includes an additional light 
field. This field is also a potential dark matter candidate with a mass m ~ 10 keV, 
that could produce another monochromatic line via its decay to photons. 

At present, there are no known models of inflation in LVS with intermediate scale 
M s . However, the Fibre Inflation model of [16| can give rise to inflation for every 



value of V. The only condition, which fixes V ~ 10 3 , and so M s ~ Mgut, is the 
matching with the COBE normalisation for the density fluctuations. Such a small 
value of V is also necessary to have a very high inflationary scale (close to the GUT 
scale) which, in turn, implies detectable gravity waves. However, in principle it is 



24 However, recently it has been discovered with the INTEGRAL spectrometer SPI [^6[ that the 511 KeV 
line emission appears to be asymmetric. This distribution of positron annihilation resembles that of low 
mass X-ray binaries, suggesting that these systems may be the dominant origin of the positrons and so 
reducing the need for more exotic explanations, such as the one presented in this paper. 
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possible that the density perturbations could be produced by another scalar field 
(not the inflaton), which is playing the role of a curvaton. In such a case, one could 
be able to get inflation also for V ~ 10 15 . In this way, both inflation and TeV scale 
SUSY would be achieved within the same model, even though gravity waves would 
not be observable. It would be interesting to investigate whether such scenarios are 
indeed realisable. 



• As derived in Section 6.2, if the volume is stabilised such that V ~ 10 15 , the decom- 
pactification temperature is rather low: 10 7 GeV. 

SV case 

In this case the volume is stabilised at smaller values of the order V ~ 10 4 , which 
allows to reproduce the standard picture of gauge coupling unification with M s ~ Mqut- 
Here are the main cosmological features of these scenarios: 

• Given that in this case V < 10 13 , all the moduli have a mass m > 10 TeV, and so 
they decay before Big-Bang nucleosynthesis. Hence these scenarios are not plagued 
by any CMP. 

• As we have already pointed out in the LV case above, smaller values of V more 



naturally give rise to inflationary models, as the one presented in [16|. Here we observe 
that the predictions for cosmological observables of Fiber Inflation were sensitive to 
the allowed reheating temperature. Since for V ~ 10 GeV we have T^ H < T max ~ 
10 15 GeV and since in |l(| the authors considered already a more stringent upper 
bound Tfjjj < 10 10 GeV (in order to avoid thermal gravitino overproduction), the 
presence of a maximal temperature does not alter the predictions of that inflationary 
scenario. 

• Fixing the volume at small values of the order V ~ 10 3 , the decompactification 
temperature turns out to be extremely high: 10 15 GeV. 

According to the discussion above, it would seem that cosmology tends to prefer smaller 
values of V. The reason is that in the SV case there is no CMP and robust models of inflation 
are known, whereas for V ~ 10 15 the light modulus suffers from the CMP and no model of 
inflation has been found yet. Interestingly enough, the lower bound on V, derived in this 
paper, suggests exactly the opposite. Namely, larger values of V are favoured since, writing 
the volume as V ~ 10 x and recalling the definition ( |6.55 ) of the parameter c, the constraint 



T* < T max rules out a relevant portion of the (x, c)-parameter space, that corresponds to 
the SV case. 

In view of this result, let us point out again that the LV case has its advantages. For 
example, the decay of the light modulus into e + e~ could contribute to explain the origin 
of the 511 KeV line. In addition, its decay to photons could produce a clean smoking-gun 
signal of LVS. Furthermore, finding a realization of inflation, that is compatible with the 
LV case, is not necessarily an unsurmountable problem. In that regard, let us note that 
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the authors of [47] proposed a model, which relates the LV to the SV case. More precisely, 
the inflaton is the volume modulus and inflation takes place at a high scale for small values 
of V. However, after inflation the modulus ends up at a VEV located at V ~ 10 15 , thus 
obtaining TeV scale SUSY. In fact, as we have already mentioned above, it could even be 
possible to realize inflation directly in the LV case. A way to achieve that would be to 



modify the the Fibre Inflation scenario of [16|, so that the density fluctuations are generated 



by a field other than the inflaton. Such curvaton-like scenarios would be very promising 
for the generation of non-gaussianities in the CMB, as well as the realization of both low 
scale inflation and low-energy SUSY. However, due to the low inflationary scale, in these 
models gravity waves will be unobservable. 

Now, even if inflation turns out not to be a problem for the LV case, there is still the 
CMP due to the presence of the light volume modulus. The results of this paper pose a 
challenge for the solution of this problem. Indeed, as we have shown in Subsection |6.3| , the 
CMP cannot be solved by diluting the volume modulus via the entropy increase caused by 
the decay of the small moduli. The reason is that the latter moduli decay before they can 
begin to dominate the energy density of the Universe. So let us now discuss in more detail 
the prospects of the other main possible solution of the CMP in LVS, namely thermal 
inflation. 

Thermal Inflation 

Thermal inflation has been studied in the literature from the field theoretic point of 



view [45]. The basic idea is that a field eft, whose VEV is much larger than its mass (and so 
is called flaton) can be trapped by thermal corrections at a false vacuum in the origin. At 
a certain temperature, its vacuum energy density can start dominating over the radiation 
one, thus leading to a short period of inflation. This period ends when the temperature 
drops enough to destabilise the local minimum the flaton was trapped in. 

Since the flaton <f> has to have a VEV {(ft) 3> m^, it is assumed that the quartic piece 
in its potential is absent. However in this way, the 1-loop thermal corrections cannot trap 
the flaton in the origin because they go like 

Vt ~ T*ml = T 2 ^, (7.1) 

and there is no quartic term in V that would give rise to a term like T 2 (ft 2 . Hence, it 
is usually assumed that there is an interaction of the flaton with a very massive field, 
say a scalar ifi, of the form gift 2 (ft 2 , where g ~ 1 so that tp thermalises at a relatively low 
temperature. At this point, a 1-loop thermal correction due to tp would give the required 
term 

V T ~ T 2 m\ = gT 2 4> 2 . (7.2) 

When eft gets a nonzero VEV, the interaction term gift 2 (ft 2 generates a mass term for ift of 
the order ~ {cj)). Hence, when (ft is trapped in the origin at high T, ift becomes very 
light. Close to the origin, the potential looks like: 

V = Vo + (gT 2 -m 2 )^ 2 + ..., (7.3) 
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where Vq is the height of the potential in the origin. A period of thermal inflation takes 
place in the temperature window T c < T < Tj n , where Tj n ~ is the temperature 

at which the flaton starts to dominate the energy density of the Universe (beating the 
radiation energy density p r ~ T 4 ) and T c ~ ms/g is the critical temperature at which the 
flaton undergoes a phase transition rolling towards the T = minimum. The number of 
e-foldings of thermal inflation is given by: 

Let us see how the above picture relates to the LVS. In the case of V ~ 10 15 , the 
modulus t s has the right mass scale and VEV to produce N e ~ 10 e-foldings of inflation, 
which would solve the CMP without affecting the density perturbations generated during 
ordinary high-energy inflation. However, in Subsection ^O] we derived the relevant 1-loop 
temperature corrections to the scalar potential and showed that they are always subleading 
with respect to the T = potential, for temperatures below the Kaluza-Klein scale. Hence, 
since thermal inflation requires the presence of new minima at finite-temperature, we would 
be tempted to conclude that it does not take place in the LVS. In fact, this was to be 
expected also for the following reason. According to the field theoretic arguments above, 
in order for thermal inflation to occur, it is crucial that the flaton be coupled to a very 
massive field ip. However, in our model there is no particle, which is heavier than the flaton 
candidate r s . It is not so surprising, then, that we are not finding thermal inflation. 

Let us now discuss possible extensions of our model that could, perhaps, allow for 
thermal inflation to occur, as well as the various questions that they raise. 

1. Since in our case r s is the candidate flaton field, the necessary if) field would have a 
mass of the order ~ {t s )M s , and so it is likely to be a stringy mode. In such 
a case, it is not a priori clear how to compute thermal corrections to Vr due to the 
presence of ip in the thermal bath. 

2. Even if we can compute Vr, it is not clear why these corrections should trap r s at 
the origin. Note, however, that this is not implausible, as the origin is a special point 
in moduli space, where new states may become massless or the local symmetry may 
get enhanced. Any such effect might turn out to play an important role. 

3. Even assuming that Vt does trap t s in the origin, one runs into another problem. 
Namely, the corresponding small cycle shrinks below M s and so we cannot trust our 
low-energy EFT. For a full description, we should go to the EFT that applies close 
to the origin. The best known examples of these are EFTs for blow-up fields at the 
actual orbifold point. In addition, one should verify that V stays constant when the 
t s cycle shrinks to zero size. 

4. When t s goes to zero, the field ip should become massless, according to the comparison 
with the field theoretic argument (if this comparison is valid). So possible candidates 
for the role of the ip field could be winding strings or Dl-branes wrapping a 1-cycle 
of the collapsing 4-cycle. 



- 53 - 



5. If tp corresponds to a winding string, the interaction of the flaton t s with tp cannot 
be seen in the EFT and it would be very difficult to have a detailed treatment of this 
issue. 

6. The field could also be a right handed neutrino, or sneutrino, heavier than r s . The 
crucial question would still be if it would be possible to see tp in our EFT description. 
In addition, one would need to write down as a function of r s and V. It goes 
without saying that this issue is highly dependent on the particular mechanism for 
the generation of neutrino masses. 

7. Besides the small modulus r s , another possible flaton candidate could be a localised 
matter field such as an open string mode. However we notice that the main con- 
tribution to the scalar potential of this field should come from D-terms, and that a 
D-term potential usually gives rise to a mass of the same order of the VEV. Hence 
it may be difficult to find an open string mode with the typical behaviour of a flaton 
field. 

In general, all of the above open questions are rather difficult to address. This poses a 
significant challenge for the derivation of thermal inflation in LVS and the corresponding 
solution of the CMP. However, let us note that the CMP could also be solved by finding 
different models of low-energy inflation, which do not rely on thermal effects. 

8. Conclusions 

In this paper, we studied how finite-temperature corrections affect the T = effective 
potential of type IIB LVS and what are the subsequent cosmological implications in this 
context. 

We showed that the small moduli and modulini can reach thermal equilibrium with the 
MSSM particles. Despite that, we were able to prove that their thermal contribution to the 
effective potential is always subleading compared to the T = potential, for temperatures 
below the Kaluza-Klein scale. As a result, the leading temperature-dependent part of the 
effective potential is due only to the MSSM thermal bath and it turns out to have runaway 
behaviour at high T. We derived the decompactification temperature T max , above which 
the T = minimum is completely erased and the volume of the internal space starts run- 
ning towards infinity. Clearly, in this class of IIB compactifications the temperature T max 
represents the maximal allowed temperature in the early Universe. Hence, in particular, 
it gives an upper bound on the initial reheating temperature after inflation: T^ H < T max . 
25 The temperature T^ H is highly dependent on the details of the concrete inflationary 
model and re-heating process, and so in principle its determination is beyond the scope of 
our paper. Nevertheless, we can compute the temperature of the Universe after the small 
moduli decay. They are rather short-lived and their decay can either be the main source of 

25 Note, however, that it may be possible to relax this constraint to a certain degree by studying the 
dynamical evolution of the moduli in presence of finite temperature corrections, in the vein of the consid- 
erations of ^] for the KKLT set-up. 



-54- 



initial reheating (in which case the temperature after their decay is exactly T^ H ) or it can 
occur during a radiation dominated epoch, after initial reheating has already taken place. 
In both cases, the resulting temperature of the Universe T* has to satisfy T* < T max , 
which implies a lower bound on the allowed values of V. We were able to derive this bound 
and show that it rules out a large range of smaller V values (which lead to standard GUT 
theories), while favouring greater values of V (that lead to TeV scale SUSY). Note though, 
that the condition T* < T max is both necessary and sufficient in the case the decay of the 
small moduli is the origin of initial reheating, whereas it is just necessary but not sufficient 
in the case the small moduli decay below T^ H . 

Finally, we discussed possible cosmological applications of our work. In particular, we 
argued that, to realize thermal inflation in this type of compactifications, one needs to go 
beyond the current effective field theory description of the closed string moduli sector. 
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A. Moduli couplings 

We shall now assume that the MSSM is built via magnetised D7 branes wrapping an 
internal 4-cycle within the framework of 4D M = 1 supergravity. The full Lagrangian of 
the system can be obtained by expanding the superpotential W, the Kahler potential K 
and the gauge kinetic functions fi as a power series in the matter fields: 

W = W mod {<p) + K<f)H u H d + XMJ^flc^C" + (A.l) 

b 

K = K^fatp) + [Z(<p,<p)H u H d + h.c] + (A.2) 

fi = ^f^ + h i (F)S. (A.3) 

In the previous expressions, if denotes globally all the moduli fields, and W modi and K mod 
are the superpotential and the Kahler potential for the moduli, which we have discussed 
in depth in Section |2[ H u and H d are the two Higgs doublets of the MSSM, and the C's 
denote collectively all the matter fields. In the expression for the gauge kinetic function 



( |A.3|) , Tmssm is the modulus related to the 4-cycle wrapped by the MSSM D7 branes, 
and hi (F) are 1-loop topological functions of the world- volume fluxes F on different branes 
(the index i runs over the three MSSM gauge group factors). Finally the moduli scaling of 
the Kahler potential for matter fields Kfj(ip,(p) and Z(ip,if), for LVS with the small cycle 
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t s supporting the MSSM, has been derived in [38] and looks like: 26 

1/3 1/3 

%(^) ~ — %(E7) and Z^,^)--?-^). (A.4) 
A.l Moduli couplings to ordinary particles 

We now review the derivation of the moduli couplings to gauge bosons, matter particles 
and Higgs fields for high temperatures T > Mew- in this case all the gauge bosons and 
matter fermions are massless. 

• Couplings to Gauge Bosons 

The coupling of the gauge bosons X to the moduli arise from the moduli dependence of 
the gauge kinetic function (|A.3[ ). We shall assume that the MSSM D7 branes are wrapping 
the small cycle 27 , and so we identify Tmssm = T a . We also recall that the gauge couplings 
of the different MSSM gauge groups are given by the real part of the gauge kinetic function, 
and that one obtains different values by turning on different fluxes. Thus the coupling of t s 
with the gauge bosons is the same for U(l), SU(2) and 577(3). We now focus on the U(l) 
factor without loss of generality. The kinetic terms read (neglecting the r s -independent 
1-loop contribution) 

We then expand t s around its minimum and go to the canonically normalised field strength 
G uu defined as 

GW = (A.6) 

and obtain 



{--gauge — G/ii/G^ , . G^uG^ . (A. 7) 



5Tj_ 

M p {t : _ 

Now by means of (|4.8j ) we end up with the following dimensionful couplings 




M P 



XG^G^, (A.8) 
$G Mi/ G^. (A.9) 



Couplings to matter fermions 



The terms of the supergravity Lagrangian which are relevant to compute the order of 
magnitude of the moduli couplings to an ordinary matter fermion ip are its kinetic and 



mass terms 28 : 



C = K^pi-fd^ + e x/2 Aff#, (A.10) 



26 Note that, in the case of more than one small cycle supporting the MSSM, these expressions would be 
more complicated. 

27 The large cycle would yield an unrealistically small gauge coupling: g 2 ~ {i~b}~' 1 ~ 1CP 10 . 
28 Instead of the usual 2-component spinorial notation, we are using here the more convenient 4-component 
spinorial notation. 
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where H is the corresponding Higgs field (either H u or H^). The moduli scaling of 

is given in ( |A,4| ), whereas e K ^ 2 = V -1 . Expanding the moduli and the Higgs around their 

VEVs, we obtain 

(A.ll) 

We now canonically normalise the ip kinetic terms {ip — ► ip c ) and rearrange the previous 
expression as 

1 5r, 1 5ti, \ 

+ m V>VWc + £-SH, (A.12) 



3(r s ) 2(r 6 )y 



where 



™* = {rs )T { i)V* > and ^ = ( (^T^yva ) ^^-( 2(^)3/2^)1/3 ) W^c- 

(A.13) 

The second term of ( [A. 12 ) does not contribute to the moduli interactions since Feynman 
amplitudes vanish for on-shell final states satisfying the equations of motion. Writing 
everything in terms of and we end up with the following dimensionless couplings 



~ ( ! ^ J *&Vo- (A.15) 

Moreover the first term in the Higgs Lagrangian ( A.13[ ) gives rise to the usual Higgs- 
fermion-fermion coupling, whereas the second term yields a modulus-Higgs-fermion-fermion 
vertex: 

^SH^c ~ f^T/s) SH ^c, (A.16) 



1 

M P V 5 /6 



^ ~ ( T7-T^ ) (A.18) 



We notice that for T > Mew the fermions are massless since (H) = 0, and so the two 
direct moduli couplings to ordinary matter particles ( [A. 14 ) and ( A.15| ) are absent. 

• Couplings to Higgs Fields 

The form of the un-normalised kinetic and mass terms for the Higgs from the super- 
gravity Lagrangian, reads: 

Cmggs = K RH d^Hd»H - K BH {ft 2 + ml) HH, (A.19) 
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where H denotes a Higgs field (either H u or H d ), and £i and m,Q are the canonically 
normalised supersymmetric /i-term and SUSY breaking scalar mass respectively. Their 
volume dependence, in the dilute flux limit, is 3£]: 

M P 



|//| ~ rriQ 



(A.20) 



VlnV 

In addition to ( |A,19|) , there is also a mixing term of the form 

Cm ggs mix = Z (8^H d d^H u + d^H d d^H u ) - KjjjjBjl (H d H u + H d H u ) , (A.21) 

with 

Bfi ~ ml. (A.22) 

However given that we are interested only in the leading order volume scaling of the Higgs 
coupling to the moduli, we can neglect the 0(1) mixing of the up and down components, 
and focus on the simple Lagrangian: 



'mi* 



M 2 -\ 
'"' 1 ' (VlnV) 2 J 

n I □ + . M ^.o \H + H\r f 



1 



K 



HH 



(VlnV)' 



M 2 P 



(VlnV)' 



H 



(A.23) 



where we have integrated by parts. We now expand K^ H and (VlnV) 2 and get: 



Ch 



1 



iggs 



-K 1 + - 



1 6t„ 



3 (t.) (n) 

II [ □ f m 2 H (l - 3^ 



H[U + m 



H 



'(n) 



H+ 



H 



(A.24) 



where Kq = (r s ) 1 / 3 (V) 2 / 3 and the Higgs mass is given by 

M P 



(A.25) 



f K()H, we obtain 



(V)ln(V)' 

Now canonically normalising the scalar kinetic terms H — * H c 
Cmggs = --[H c (D + m 2 H ) H c + H C (D + m 2 H ) H c ] (A.26) 

The second term in the previous expression does not contribute to scattering amplitudes 
since Feynman amplitudes vanish for final states satisfying the equations of motion. Thus 
the dimensionful moduli couplings to Higgs fields arise only from the third term once we 
express 5rb in terms of <£ and \ using ([4.7|). The final result is 



m 



H 



MpVV 



§H c Hc 



( M P 



C 



XH C H C 



m 



H 



M P 



XH C H C 



\V 5 /2(lnV) 
M P 



<5>H C H C , 



V 2 (lnV) 2 



XH C H C . 



(A.27) 
(A.28) 
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A. 2 Moduli couplings to supersymmetric particles 

We shall now work out the moduli couplings to gauginos, SUSY scalars and Higgsinos. 
Given that we are interested in thermal corrections at high temperatures, we shall focus 
on T > Mew- Thus we can neglect the mixing of Higgsinos with gauginos into charginos 
and neutralinos, which takes place at lower energies due to EW symmetry breaking. 

• Couplings to Gauginos 



The relevant part of the supergravity Lagrangian involving the gaugino kinetic terms 
and their soft masses looks like 



Cgaugino ^ -^-X'ia^X' + ^ (A' A' + h.C.) , 



(A.29) 



where in the limit of dilute world-volume fluxes on the D7-brane, the gaugino mass is given 
by M1/2 = Jf~ (Hi- Now if the small modulus supporting the MSSM is stabilised via 
non-perturbative corrections, then the corresponding F-term scales as 



F * * vEv" (A ' 30) 

Notice that the suppression factor InV ~ ln(Mp/m 3/ / 2 ) in (|A.30|) would be absent in the 
case of perturbative stabilisation of the MSSM cycle Q. Let us expand t s around its VEV 
and get: 



'"gaugino — 



(Ts) 



X'ia^X' + \^(X'X' + h.c.) 



+ 



5ts_ 
M P 



X'ia^daX' + 



M P (A' A' + h.c. 



11 ' ' (V) ln(V) 2 

(A.31) 

We need now to expand also r& around its VEV in the first term of ( A.31] ): 



VlnV 



3/2 



lnV (V)ln(V)V 2<r 6 ) 



1--^ + ... 



(A.32) 



and canonically normalise the gaugino kinetic terms A' — > A = \J (r s )A'. At the end we 
obtain: 



■jnugmo — AfO Up AT 



M P {XX + h.c.) {XX + h.c.) ( 5t s 3 8ti 



+ 2<V)ln(V> Ws) 2 {n) ) + (r s )Mp A ^ A ' 



From ([AT33D we can immediately read off the gaugino mass: 



Mi 



M F 



F s 



rn 



3/2 



1/2 (V) ln(V) r s ~ In (Mp/m 3/2 ) 
Let us now rewrite ( |A.33 ) as: 



(A.33) 
(A.34) 



C 



gaugino 



- I 1 + ) 



Mi/a 

Xia^daX + {XX + fc.c. 



3 s n M 1/2 
4(r 6 ) M P 



(AA + /».&). 



(A.35) 



- 59 - 



We shall now focus only on the last term in ( |A.35| ) since it is the only one that contributes 
to decay rates. In fact, Feynman amplitudes with on-shell final states that satisfy the equa- 
tions of motion, are vanishing. Using Q4.7D , we finally obtain the following dimensionless 
couplings: 



£*AA 



M- 



1/2 



MpVV 



$AA 



1 



£ 



M 



1/2 



M P 



XAA 



V 3 / 2 InV 
1 



<3?AA, 



VlnV 



XAA. 



(A.36) 
(A.37) 



• Couplings to SUSY Scalars 

The form of the un-normalised kinetic and soft mass terms for SUSY scalars from the 
supergravity Lagrangian, reads: 



£ 



scalars 



K 



a 3 



(VlnV) z 



Assuming diagonal Kahler metric for matter fields 



— Ka5 al3 , 



(A.38) 



(A.39) 



the initial Lagrangian ( A.38| ) simplifies to 



£ scalars — K a yd^C a d^C C 
= ~Ka C" \ □ 



(VlnV) 
1 



C a c c 



c a + c a □ 



(A.40) 



(VlnV) 2 / ' V (VlnV) 2 , 

We note that ( [A. 40 ) is of exactly the same form as the Higgs Lagrangian (A. 23). This 
is not surprising since for temperatures T > M^W: tli e Higgs behaves effectively as a 
SUSY scalar with mass of the order the scalar soft mass: tuh ~ rriQ. Thus we can read 
off immediately the dimensionful moduli couplings to the canonically normalised SUSY 
scalars tp from ( |A.27| ) and ( |A.28| ): 



M P 



( M P 



„ ( rrin 

Imp ' V " 



V 2 (lnV) 2 



(A.41) 
(A.42) 



• Couplings to Higgsinos 

The relevant part of the supergravity Lagrangian involving the Higgsino kinetic terms 
and their supersymmetric masses looks like: 



Ch 



iggsmo — HH 



Huia^d^Hu + H d ia^d^H d + p, H u H d + h.c 



(A.43) 
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After diagonalising the supersymmetric Higgsino mass term, we end up with a usual La- 
grangian of the form: 



iggsmo — 1Y HH 



Hia^duH + y.[HH + h.c. 



(A.44) 



where H denotes collectively both the Higgsino mass eigenstates, which are the result of 
a mixing between the up and down gauge eigenstates. We recall also that since we are 
focusing on temperatures above the EWSB scale, we do not have to deal with any mixing 
between Higgsinos and gauginos to give neutralinos and charginos. Expanding the Kahler 
metric (A4) and the /U-term ( A. 20 ), we obtain: 

1 5t s 5r b 



Ch 



iggsmo — ^0 



K 1 + 



3 (t s ) (n) 



HH + h.c. 



(A.45) 

where Kq = (t s ) 1//3 (V)~ 2 / 3 and the physical Higgsino mass is of the same order of magni- 
tude of the soft SUSY masses: 



M P 



Now canonically normalising the scalar kinetic terms H — > H, 

1 6t s Sn \ 



(A.46) 

'KqH , we end up with: 



iggsmo 



1 + 3(r s ) (r b )J 

3 5n 



H c ia»d a H c + (H C H C + h.c. 



(A.47) 



Writing everything in terms of $ and Xi from the last term of ( A.47 ), we obtain the 
following dimensionless couplings: 



xHcH c 



in 



H 



M P 



in 



XH C H C 



1 



VlnV 



XH C H C 



MpVV 



$H C H C 



1 



V 3 / 2 lnV' 



(A.48) 
(A.49) 



A. 3 Moduli self couplings 

In this Section we shall investigate if moduli reach thermal equilibrium among themselves. 
In order to understand this issue, we need to compute the moduli self interactions, which 
can be obtained by first expanding the moduli fields around their VEV 

T i = (T i ) + dT i , (A.50) 

and then by expanding the potential around the LARGE Volume vacuum as follows: 



1 d 2 V 

v = v((r s ),(T b )) + - 



2 dTidrj 



1 



d 3 V 



i^Tj + ^ dTidTjChfc 



5Ti5Tj5T k + .. 



(A.51) 



We then concentrate on the trilinear terms which can be read off from the third term of 
( A.51| ). We neglect the 0{5rf) terms since the strength of their couplings will be subleading 
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with respect to the 0{5rf) terms since one has to take a further derivative which produces a 
suppression factor. Taking the third derivatives and then expressing these self-interactions 
in terms of the canonically normalised fields 

5n ~ o (v 1/6 ) $ + o (v 2/3 ) x , 

we end up with the following Lagrangian terms at leading order in a large volume expansion: 

£*3 c ^\ C^ x c x * 2 , (A.52) 

Mp n Mp n 

C x^ = > ^ x a ~ ^ x 3 . (A.53) 
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